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ABSTRACT 

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The 
catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy 
clusters. It contains 1653 detections, of -which 1203 are confirmed clusters with identified counterparts in external data-sets, and 
is the first SZ-selected cluster survey containing > 10® confirmed clusters. We present a detailed analysis of the survey selection 
function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find 
that the Ysrsoo estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires, 
the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, 
which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the 
counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- 
luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their 
mass, but are almost absent from ROSAT X-ray selected samples. 

Key words, cosmology: observations - galaxies: clusters: general - catalogues 
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1. Introduction 

This paper is one of a set associated with the 2015 Planck^ 
full mission data release and describes the production and 
properties of the legacy catalogue of Sunyaev Zeldovich 
sources (PSZ2). 

In the framework of hierarchical structure formation, 
peaks in the cosmological density field collapse and merge 
to form gravitationally bound haloes of increasing mass 
(Peebles 1980). The galaxy clusters are the most massive of 
these bound structures and act as signposts for the extrema 
of the cosmological density field on the relevant scales. The 
evolution of galaxy cluster abundance with mass and red- 
shift is thus a sensitive cosmological probe of the late-time 
universe, providing unique constraints on the normalisa¬ 
tion of the matter density fluctuations, (Tg, the mean mat¬ 
ter density, flm, the density and equation of state of the 
dark energy field, Ude and w, as well as constraining some 
extensions of the minimal cosmological model, such as mas¬ 
sive neutrinos, and non-standard scenarios such as modified 
gravity (see eg: Borgani & Kravtsov 2009; Allen et al. 2011). 
In recent years, cluster data from the microwave through to 
the X-ray have been used to constrain cosmological param¬ 
eters (Vikhlinin et al. 2009b; Rozo et al. 2010; Hasselfield 
et al. 2013; Benson et al. 2013; Planck Collaboration XX 
2014; Zu et al. 2014). 

Galaxy clusters are multi-component objects composed 
of dark matter, which dominates the mass, stars, cold 
gas and dust in galaxies, and a hot ionised intra-cluster 
medium (ICM). These different components make clusters 
true multi-wavelength objects. The galaxies emit in the op¬ 
tical and infrared. The ICM, which is the majority of the 
baryonic material by mass, emits in the X-rays via thermal 
bremsstrahlung and line emission, and energy-boosts cos¬ 
mic microwave background photons via inverse Compton 
scattering. 

This last effect, the thermal Sunyaev Zeldovich (SZ) ef¬ 
fect (Sunyaev & Zeldovich 1970, 1980), imprints a redshift- 
independent spectral distortion on the cosmic-microwave 
background (CMB) photons reaching us along the line of 
sight to the cluster. This results in an increase in intensity 
at frequencies above 220 GHz, and a decrease in intensity at 
lower frequencies. The Planck High-Frequency Instrument 
(HFI) is unique in providing high-precision data for both 
the increment and the decrement across the whole sky. 

The utility of a cluster survey for cosmological work de¬ 
pends on our ability to determine accurately its selection 
function and to obtain unbiased measurements of cluster 
mass and redshift. The first cluster surveys consisted of 
galaxy overdensities identified by eye from photographic 
plates (Abell 1958). The construction of large optical cata¬ 
logues improved significantly with the data from the SDSS 
(Koester et al. 2007), whose five photometric bands have 
allowed robust photometric classification of red-sequence 
cluster galaxies and accurate photometric redshifts to z < 
0.6 across 1/4 of the sky (Hao et al. 2010; Szabo et al. 
2011; Wen et al. 2009; Rykoff et al. 2014). These catalogues 

^ Planck (http://www.esa.int/Plaiick) is a project of the 
European Space Agency (ESA) with instruments provided by 
two scientific consortia funded by ESA member states and led 
by Principal Investigators from France and Italy, telescope re¬ 
flectors provided through a collaboration between ESA and a 
scientific consortium led and funded by Denmark, and additional 
contributions from NASA (USA). 


now typically contain 10 "^ — 10 ® clusters and provide cluster 
richness as an observable that correlates with mass with an 
intrinsic scatter (Tint of about 25% (Rozo & Rykoff 2014). 

Construction of X-ray cluster surveys is now a mature 
activity, with several catalogues now available based on all¬ 
sky data from the ROSAT satellite, alongside additional 
catalogues of serendipitous detections from pointed obser¬ 
vations (Ebeling et al. 1998; Bohringer et al. 2004; Reiprich 
& Bohringer 2002; Ebeling et al. 2010; Piffaretti et al. 2011; 
Burenin et al. 2007; Mehrtens et al. 2012). The most basic 
X-ray survey observable, the X-ray luminosity L 500 mea¬ 
sured within rsoo has been shown correlate with mass 
with intrinsic scatter of about 40% (Pratt et al. 2009). 
Observables with lower intrinsic scatter against mass can be 
dehned when pointed X-ray follow-up information is avail¬ 
able, including the core-excised X-ray luminosity (Maughan 
2007; Pratt et al. 2009) and Fx, the product of the gas mass 
and the core-excised spectroscopic temperature (Kravtsov 
et al. 2006; Vikhlinin et al. 2009a; Mahdavi et al. 2013). 
While X-ray surveys are unique in their purity, they do 
suffer from selection biases that favour low-redshift sys¬ 
tems, due to flux limitations, and dynamically relaxed clus¬ 
ters with an X-ray bright cooling core (Eckert et al. 2011; 
Schuecker et al. 2003; Vikhlinin et al. 2009b; Chen et al. 
2007). 

SZ surveys offer a different window on the cluster popu¬ 
lation: their selection function flattens towards higher red- 
shifts, providing a nearly mass-limited census of the cluster 
population at high redshift, where abundance is strongly 
sensitive to cosmological parameters (Carlstrom et al. 2002; 
Planck Collaboration XXIX 2014). The SZ survey observ¬ 
able is the spherically integrated Comptonisation parame¬ 
ter, Ysz , which is related to the integrated electron pressure 
and hence the total thermal energy of the cluster gas. It 
is also expected to correlate with mass with a low intrinsic 
scatter and little dependence on the dynamical state of the 
cluster (eg: da Silva et al. 2004; Kay et al. 2012; Hoekstra 
et al. 2012; Planck Collaboration Int. HI 2013; Sifon et al. 
2013). 

The spherically-integrated pressure profiles of X-ray and 
SZ clusters have been observed to follow a near universal 
profile with little dispersion (Arnaud et al. 2010; Planck 
Collaboration Int. V 2013), permitting the detection of 
clusters with a matched multi-frequency filter based on 
some assumed pressure profile (Herranz et al. 2002; Melin 
et al. 2006a). Samples constructed this way have well un¬ 
derstood selection functions, though discrepancies due to 
profile mismatch or contaminating infra-red emission may 
still be present to some level. Large SZ surveys have only 
appeared recently, with catalogues of order ~ 10 ^ clusters 
released by the Atacama Cosmology Telescope (Hasselfield 
et al. 2013), the South Pole Telescope (Reichardt et al. 
2013) and Planck satellite collaborations. 

This is the third all-sky catalogue produced from Planck 
SZ data. The early Sunyaev-Zeldovich (ESZ) catalogue pre¬ 
sented 189 clusters detected from 10 months of survey data 
(Planck Collaboration 2011), while the PSZl, the full-sky 
catalogue assembled from the nominal mission data, pre¬ 
sented 1227 cluster candidates detected from 15.5 months 
of data (Planck Collaboration XXIX 2014). This paper 


^ T 500 is the cluster-centric distance within which the mean 
density is 500 times the critical density of the Universe at the 
cluster redshift. 
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presents 1653 candidates detected from the full mission 
survey of 29 months. 1203 of these have been confirmed 
in ancillary data, and 1094 have redshift estimates. The 
PSZ2 expands the scope and sensitivity of the SZ view of 
galaxy clusters by substantially increasing the number of 
lower mass clusters available for study. It is also expected to 
contain many new, as yet unconfirmed, high-redshift clus¬ 
ters. We report on the construction and characterisation 
of the catalogue, presenting the survey selection functions 
and a compilation of multi-wavelength ancillary informa¬ 
tion including redshifts. We also briefly discuss the physical 
properties of the sample. 

This paper is organised as follows. In Sect. 2 we sum¬ 
marise the three extraction algorithms used to build the 
catalogue, focussing on the changes in the algorithms since 
they were used to construct the PSZl. In Sect. 3 we describe 
the construction of the catalogue. In Sect. 4 we present the 
survey selection functions (completeness and statistical re¬ 
liability) and the complementary approaches used to esti¬ 
mate them, while in Sect. 5 we discuss and validate the 
estimation of the Ysz parameters, both blindly and when 
using prior information, and we compare the consistency of 
the new catalogue with the PSZl in Sect. 6. In Sect. 7 we 
report on the search for multi-wavelength counterparts in 
ancillary catalogues and follow-up observations. Finally we 
present the physical properties of the sample in Sect. 8 and 
conclude in Sect. 9. A full description of the available data 
products is given in Appendix D 

2. Extraction Algorithms 

The SZ detection and parameter estimation algorithms 
used to construct the PSZ2 extend and refine those used to 
construct the PSZl. In this section we recall the principles 
of the three algorithms. The refinements of each algorithm 
since the PSZl release are detailed in Appendix C. Two of 
the algorithms (MMFl and MMF3) are based on the same tech¬ 
nique (Matched Multi-filters) but have been implemented 
independently^. The third one (PwS for PowellSnakes) relies 
on Bayesian inference. 

2.1. Matched Multi-filters: MMFl and MMF3 

The matched filtering technique was first proposed for SZ 
studies by Haehnelt & Tegmark (1996). It was subsequently 
developed by Herranz et al. (2002) and Melin et al. (2006b) 
for SZ cluster extraction in multifrequency data sets such 
as Planck. The method was later adopted by the SPT 
and ACT collaborations (Staniszewski et al. 2009; Marriage 
et al. 2011). 

We model the vector of map emission at each frequency 
m{x), at a given position on the sky x as 

m{x) = yote, (x) -F n{x) (1) 

where tg, {^) is the signal vector describing the spatial dis¬ 
tribution at each frequency of the SZ emission from a clus¬ 
ter with angular size 9s, n{x) is the total astrophysical and 

^ The MMF numbers were given after the comparison of 
twelve algorithms in an earlier phase of the Planck mis¬ 
sion (Melin et al. 2012). MMFl and MMF3 were respectively the 
first and third algorithm based on Matched Multi Filters to en¬ 
ter the comparison 


instrumental noise. The frequency component of the sig¬ 
nal vector is the normalized cluster profile (a;) (Arnaud 
et al. 2010) convolved by the Planck beams bi{x) and scaled 
with the characteristic frequency dependance of the 

thermal SZ effect: tg^{x)i = jy{i'i)\bi*Te^\{x). 9s is the clus¬ 
ter scale radius, which is related to 05 oo through the con¬ 
centration parameter C 500 by 6*500 = C 500 x 9s- The Matched 
Multi-filter allows us to recover an unbiased estimate 
ijo of the central Comptonization parameters j/o with min¬ 
imal variance : 

yo = J d^x'i>gj{x)-m{x), ( 2 ) 

where 

^eAk)=<TlP-\k)-teAk), (3) 

with 

d'^k teJik) ■ P~^ ■ teAk) 

P{k) being the cross-channel power spectrum matrix of 
the maps. It is effectively the noise matrix for the MMF, 
because the tSZ is small compared to other astrophysical 
signals, and is estimated directly from the maps. 

The MMF algorithms first divide each Planck all-sky 
map in 640/504 tangential maps (14.66/10 degrees on a 
side) for MMFl /MMF3 respectively. Each set of tangential 
maps is Altered by with the assumed cluster size vary¬ 
ing from 6*;,=0.8 to 32 arcmin. We then locate peaks in the 
filtered maps above a S/N threshold of four. The locations 
of the peaks give the positions of our cluster candidates. 
These are then merged into a single all-sky catalogue by 
merging candidates separated by less than 10 arcmin. For 
MMF3, we performed a second step by creating sets of smaller 
rectangular frequency maps centered on the cluster candi¬ 
dates identified in the first step. We re-apply the MMF on 
these centred tangential maps which allows a better estima¬ 
tion of the background. During the second step, the sizes 
and fluxes are estimated more precisely. This second step 
is only performed for MMF3 because the overlap of the tan¬ 
gential maps in the first step is small with compared to 
MMF land PwS. 

We define the blind cluster size as the filter scale that 
maximizes the S/N at the location of the cluster candidate 
and the blind flux is defined as the corresponding 1/0 pa¬ 
rameter. We then define the integrated blind flux as: 

Y 5 R 500 = 2/0 / drrg^ (r) (5) 

J ^<5x^500 

Each of the algorithms produces probability distribu¬ 
tions in the {9s, F 5 R 500 ) plane for each detection, marginal¬ 
ising over the parameters for the centre of the cluster, which 
possess a Gaussian likelihood. The algorithms also return 
an estimate of the radial position uncertainty, derr from the 
position likelihood. 

Although the two implementations of the MMF are 
quite close, they produce noticeably different catalogues be¬ 
cause the extraction is very sensitive to the estimation of 
the background (Eq. 4). Both the size adopted for the tan¬ 
gential maps and the details of the estimation of the matrix 
P{k) impact the S/N and hence which peaks are detected. 
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Table 1. Effective frequencies and Gaussian beam widths 
assumed for extraction per channel. The beam widths 
are the mean full-width-at-half-maximum of the FEBeCoP 
Gaussian beam fits across the sky, in arcmin. The effec¬ 
tive frequencies t'eff, shown in GHz, encapsulate band-pass 
effects in each channel. 


Channel 

FWHM 

t'eff 

100 

9.659 

103.416 

143 

7.220 

144.903 

217 

4.900 

222.598 

353 

4.916 

355.218 

545 

4.675 

528.400 

857 

4.216 

776.582 


2.2. PowellSnakes (PwS) 

PowellSnakes (PwS) is a fast, fully Bayesian, multifrequency 
detection algorithm designed to identify and characterize 
compact objects buried in a diffuse background as described 
in Carvalho et al. (2009, 2012). PwS operates using about 
2800 square patches of 14.66 degree on a side, in order to en¬ 
sure highly redundant sky coverage. PwS detects candidate 
clusters and at the same time computes the evidence ratio 
and samples from the posterior distributions of the cluster 
parameters. Then, it merges the sub-catalogues from each 
patch map and applies criteria for acceptance or rejection of 
the detection, as described in Carvalho et al. (2012). Priors 
may be provided for the position, integrated flux and radius 
of the clusters. For cluster detection, we apply flat priors 
on the position and non-informative priors in the radius 
and integrated flux, as determined using Jeffrey’s method. 
PwS uses a calibration of the cross-power spectrum that 
uses an iterative scheme to reduce the contamination of 
the background by the SZ signal itself. This makes PwS 
particularly robust to small changes in the background. 

3. Catalogue Construction 

The main catalogue is constructed by combining the detec¬ 
tions made by the three methods into a union catalogue, 
while merging the detections made by multiple more than 
one method. Half of the detections in this union are also in 
the intersection catalogue, defined as those detections made 
by all three codes simultaneously. This section describes the 
technical details of the construction of these catalogues. 

3.1. Pipeline 

The SZ catalogue construction pipeline is shown in 
schematic form in Fig. 1 and largely follows the process used 
to build the PSZl. The Planck data required for the con¬ 
struction of the catalogue comprises the HFI maps, point 
source catalogues for each of the HFI channels, effective 
frequencies and beam widths per HFI channel as shown in 
Table I, survey masks based on dust emission as seen in 
the highest Planck channels, and the catalogue of extended 
galactic cold-clump detections. 

The HFI maps are pre-processed to fill areas of miss¬ 
ing data (typically a few pixels), or areas with unusable 
data, specifically bright point sources. Point sources with 
S/N > 10 in any channel are masked out to a radius of 
ScTbeam, using a harmonic infilling algorithm. This prevents 



Fig. 1. Pipeline for catalogue construction. 


spurious detections caused by Fourier ringing in the filtered 
maps used by the detection algorithms. As a further guard 
against such spurious detections, we reject any detections 
within 5(Tbeam of a hlled point source. We have verified that 
this treatment reduces spurious detections due to bright 
point sources to negligible levels in simulations, while re¬ 
ducing the effective survey area by just 1.4% of the sky. 
Together with the 15% galactic dust and Magellanic cloud 
mask, this defines a survey area of 83.6% of the sky. 

After infilling, the three detection codes produce indi¬ 
vidual candidate catalogues down to a threshold S/N > 4.5. 
The catalogues are then merged to form a union catalogue, 
using the dust and extended point source masks discussed 
above to define the survey area. The merging procedure 
identifies the highest S/N detection as the reference posi¬ 
tion during the merge: any detections by other codes within 
5 arcmin are identified with the reference position. The ref¬ 
erence position and S/N are reported in the union cata¬ 
logue. 

PwS can produce a small number of high-significance 
spurious detections associated with galactic dusty emission. 
We apply an extra cut of PwS-only detections at S/N > 10 
where the spectrum has a poor goodness-of-fit to the SZ 
effect, > 16. 

We also remove five PSZ2 detections that match 
PSZl detections confirmed to be spurious by the 
PSZl follow-up program (these were the ones that we re¬ 
detected: there were many more confirmed spurious detec¬ 
tions from the program). 

Finally the sample is flagged to identify the various sub¬ 
samples discussed in Sect. 3.3. The most important of these 
flags is discussed in the next section. 

3.2. Infra-red spurious detections 

Gold compact infra-red emission, particularly that due to 
galactic cold-clumps, can lead to high-significance spurious 
detections. We identify these detections by searching for 
7 arcmin matches with the Planck cold-clump catalogue 
(G3PO), or with PCCS2 detections at both 545GHz and 
857GHz. This matching radius was chosen because it is 
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Fig. 2. The distribution of raw SZ detections, with deleted 
infra-red flagged candidates in red and retained infra-red 
flagged detections in green. 


the typical size of a Planck detected cold-clump (Planck 
Collaboration XXVIII 2015). 

318 raw union detections match these criteria. They 
tightly follow the distribution of galactic emission (see Fig. 
2), such that if the 65% galactic dust mask (used for clus¬ 
ter cosmology ) is used to define the sample instead of the 
85% dust mask, the number of IR-matched raw detections 
drops to 40. For the high-purity sample formed from the 
intersection of all three codes, the numbers are 82 and 13 
for the 85% and 65% dust masks respectively. Some high 
latitude spurious candidates remain. To minimise the effect 
of these probably spurious detections on the catalogue, we 
delete them. 

We have retained in the sample all 15 confirmed clusters 
that match these criteria. These IR contaminated clusters 
represent about 1.5% of the total confirmed clusters in the 
PSZ2. In the catalogue, we define a flag, IR_FLAG, to de¬ 
note the retained clusters that match these criteria. They 
can be expected to have heavily contaminated SZ signal. 

A small fraction of the unconfirmed detections deleted 
due to IR-contaminations may have been real clusters. 
Assuming that optical and X-ray confirmation is unbiased 
with regard to the presence of IR emission, we estimate 
these deletions to bias our completeness estimates by less 
than 1%. 

3.3. Catalogue sub-samples 

The union catalogue can be decomposed into separate sub¬ 
samples, defined as the primary catalogues of the three in¬ 
dividual detection codes (PwS,MMFl,MMF3), as well as into 
unions and intersections thereof. The intersection subsam¬ 
ple of candidates detected by all three algorithms can 
be used as a high-reliability catalogue with less than 2% 
spurious contamination outside of the galactic plane (see 
Sect. 4.6). 


3.3.1. The Cosmology Catalogues 

We constructed two cluster catalogues for cosmology stud¬ 
ies from the MMF3 and the intersection sub-samples respec¬ 
tively. For these catalogues, our goal was to increase as 
much as possible the number of detections while keeping 
contamination negligible. A good compromise is to set the 
S/N threshold to 6 and apply a 65% galactic and point 
source mask as in our 2013 cosmological analysis (Planck 
Collaboration XX 2014). In this earlier paper, our baseline 
MMF3 cosmological sample was constructed using a thresh¬ 
old of 7 on the 15.5 month maps, which is equivalent to 8.5 
on the full mission maps. Estimations from the QA (Fig. 11) 
suggest that our 2014 intersection sample should be > 99% 
pure for a threshold of 6. 

The MMF3 cosmological sample contains 439 detections 
with 433 confirmed redshifts. The intersection cosmologi¬ 
cal sample contains 493 detections with 479 confirmed red- 
shifts. Assuming that all detections having VALIDATION 
flag greater than zero are clusters, the empirical purity of 
our samples are > 99.8% for MMF3 and > 99.6% for the in¬ 
tersection. Note that the intersection sample contains more 
clusters than the MMF3 sample for the same S/N threshold. 
This is expected since the definition of the S/N for the in¬ 
tersection sample is to use the highest value from the three 
detection methods. 

The completeness is also a crucial piece of information. 
It is computed more easily with the single method cata¬ 
logue for which the analytical error-function (ERF) approx¬ 
imation can be used (as defined in Planck Collaboration 
XXIX 2014). In Sect. 4.3 and in Planck Collaboration XXIV 
(2015), we show that this analytical model is still valid for 
the considered threshold. For the intersection sample, we 
rely on the Monte-Carlo estimation of the completeness de¬ 
scribed in Sect. 4.2. 

These two samples are used in the cosmology analysis 
of Planck Collaboration XXIV (2015). Detections that are 
included in either of the cosmology samples are noted in 
the main catalogue (see Appendix D). 

3.4. Consistency between codes 

We construct the union sample using the code with the 
most significant detection to supply the reference position 
and S/N. This contrasts with the PSZl, which used a pre¬ 
defined code ordering to select the reference position and 
S/N. In this section, we demonstrate the consistency of the 
detection characteristics of the codes for common detec¬ 
tions, motivating this change in catalogue construction. 

We fit the S/N relation between codes using the 
Bayesian approach described by Hogg et al. (2010) for lin¬ 
ear hts with covariant errors in both variables. We consider 
the catalogue S/N values to be estimates of a true underly¬ 
ing variable, s, with Caussian uncertainties with standard 
deviation cr = 1. 

We relate the s values for two different catalogues using 
a simple linear model 

S2 = asi -F A, (6) 

where we assume flat priors for the intercept A, a flat prior 
on the arc-tangent of the gradient a, such that p{a) oc 
1/(1 -F a^). We also allow for a Caussian intrinsic scatter 
between the s values that includes any variation beyond 
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Fig. 3. Comparison of the S/N estimates from the three detection codes. The dashed green curves show the best ht 
relation for 0.8 correlation and the red line is the line of equality. 


Table 2. Results of fits between S/N from the three detection codes, using the fitting function in equation 6 . The assumed 
correlation of the uncertainties of si and S 2 was 0 . 8 . 


Si 

Si 

N 

A 

a 

(Tint 

MMFl 

MHF3 

1032 

-0.01 ±0.01 

1.01 ±0.01 

0.033 ± 0.001 

MMFl 

PwS 

985 

-0.02 ±0.01 

1.03 ±0.02 

0.030 ± 0.001 

MMF3 

PwS 

1045 

0.0 ±0.01 

1.01 ±0.01 

0.031 ±0.001 


the measurement uncertainty on s. This is parameterised 
by (Tint with an uninformative prior p{crint) oc l/crint- 

We assume a hducial correlation of Pcorr = 0.8 be¬ 
tween the S/N estimates of each code pair, which is typi¬ 
cal of the correlation between the matched-multifrequency- 
filtered patches of each code. The fit results are shown in 
Table 2. 

The S/N estimates from the three codes are compared 
in Fig. 3, which also shows the best fit relation. MMFl pro¬ 
duces noticeably weaker detection than the other two codes 
for the 14 very strong detections at S/N > 20. Excluding 
these exceptional cases from the comparison, the best-fit 
relations between the S/Ns from each code show no signif¬ 
icant deviations from equality between any of the codes. 

There are a small number of highly significant outliers 
in the relation between PwS and the MMF codes. These 
are clusters imbedded in dusty regions where the different 
recipes for the filtered patch cross-power spectrum vary sig¬ 
nificantly and the likelihood assumptions common to all 
codes break down. PwS shows outlier behaviour relative to 
the other codes as its recipe is most different from the other 
codes. 

Fig. 4 shows the consistency of the position estimates 
between the codes. The positions of MMFl and MMF3 are 
more inconsistent with one another than any other code 
combination. The 67% bound on the MMFl- MMF3 separation 
is 1.34 arcmin, while for MMFl- PwS it is 0.98 arcmin and for 
MMF3- PwS it is 1.1 arcmin. This is consistent with the obser¬ 
vation from the quality assessment that the PwS positions 
are the most robust (Sect. 4.4). 

4. Selection Function 

A necessary element of any cluster sample is the selection 
function that relates the detected sample to the underlying 
population of objects. The selection function comprises two 
complementary functions: the completeness, which defines 
the probability that a given real object will be detected; 


o 



Separation [arcmin] 

Fig. 4. Cumulative distribution of angular separation be¬ 
tween matched detections for each possible code pair. The 
vertical dashed line indicates the width of a Healpix pixel 
at the Planck resolution. 


and the statistical reliability, also known as purity, which 
dehnes the probability that a given detection corresponds to 
a real object. As a function of underlying object attributes, 
the completeness is a function of underlying SZ observables, 
6*500 and Ysoo. The reliability is a statistical function of de¬ 
tection attributes and is presented as a function of detection 
S/N. 

4.1. Monte-Carlo Injection 

The selection function is determined by the Monte-Carlo 
injection of simulated clusters into both real and simulated 
Planck maps. A common segment is the injection of cluster 
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Fig. 5. The 910 simulated pressure profiles from the cosmo- 
OWLS simulations used for cluster injection. Also shown 
are the assumed extraction profile (UPP) and the best-fit 
profile from a sample of 62 pressure profiles fitted using 
Planck and x-ray data (PIPV, Planck Collaboration Int. V 
2013). 


SZ signal. The cluster signal is assumed to be spherically 
symmetric and to follow a pressure profile similar to the 
generalised Navarro-Frenk-White (GNFW) profile assumed 
in the catalogue extraction. 

To include the effects of system-on-system variation in 
the pressure distribution, we draw the spherically-averaged 
individual pressure profiles from a set of 910 pressure pro¬ 
files from simulated clusters from the cosmo-OWLS simu¬ 
lations (Le Brun et al. 2014; McCarthy et al. 2014), an ex¬ 
tension of the Overwhelmingly Large Simulations project 
(Schaye et al. 2010). These pressure profiles are empirical 
in the sense that they have not been fitted using a GNFW 
profile: the mean pressure is used within concentric radial 
shells (after the subtraction of obvious sub-structures) and 
the injected profiles are interpolated across these shells. The 
simulated clusters were selected for this sample by requir¬ 
ing that their mass be above the approximate limiting mass 
for Planck at that redshift. The ensemble of simulated pro¬ 
files are shown in Fig. 5. Each profile is normalised such 
that the spherically integrated Y parameter matches the 
fiducial injected (^ 500 ,^ 500 ) parameters for the halo. The 
injected (^ 500 ,^ 500 ) are different for completeness and reli¬ 
ability simulations and each is discussed below. 

Effective beam variation is an important consideration 
for the unresolved clusters at the intermediate and high red- 
shifts of cosmological interest. The injected clusters are con¬ 
volved with effective beams in each pixel including asym¬ 
metry computed following Mitra et al. (2011) 

4.2. Completeness 

The completeness is defined as the probability that a cluster 
with a given set of true values for the observables (¥ 500 ,^ 500 ) 
will be detected, given a set of selection criteria. A good 
approximation to the completeness can be defined using 
the assumption of Gaussianity in the detection noise. In 
this case, the completeness for a particular detection code 
follows the error function (ERF), parametrised by a se¬ 


lection threshold q and the local detection noise at the 
clusters radial size, crY(^,'0) (see the discussion in Planck 
Gollaboration XXIX 2014). This approach is not suited to 
the union and intersection catalogues from Planck due to 
the difficulty in modelling correlations between detection 
codes. We determine the completeness by brute force: in¬ 
jecting and detecting simulated clusters into the Planck sky 
maps. This approach has the advantage that all algorithmic 
effects are encoded into the completeness, and the effects of 
systematic errors such as beam and pressure profile varia¬ 
tion can be characterised. This approach also fully accounts 
for the non-Gaussianity of the detection noise due to fore¬ 
ground emission. 

The injected (Vsooj^soo) parameters are drawn from a 
uniform distribution in the logarithm of each variable, en¬ 
suring that our logarithmically spaced completeness bins 
have approximately equal numbers of injected sources. 

As the completeness is estimated from injection into real 
data, injected sources can contribute to the detection noise. 
We therefore use an injection mode, as was the case for 
the PSZl completeness, where injected clusters are removed 
from the maps used to estimate the noise statistics. We 
also avoid superimposing injected clusters on top of one 
another, or on top of real data detections. Together, these 
ensure that the noise statistics for injected clusters are the 
same as for the real detections in the map. 

We release as a product the Monte-Carlo completeness 
of the catalogues at thresholds stepped by 0.5 in S/N over 
the range 4.5 <S/N< 10. Fig. 6 shows the completeness of 
the union and intersection catalogues as functions of input 
(^sooifi'soo) and at representative values of O 500 , for three 
detection thresholds. The union and intersection catalogues 
are most similar at high S/N, where they match well except 
at small scales. Here the intersection catalogue follows the 
lower completeness of MMFl. This is due to an extra selection 
step in that code which removes spurious detections caused 
by point sources. The union and the intersection catalogues 
mark the upper and lower limits of the completeness values 
for the sub-catalogues based on the individual codes. 

The completeness of the Planck cluster catalogue is ro¬ 
bust with respect to deviations of the real SZ profiles of 
galaxy clusters from the one assumed by the algorithms for 
filter construction. To demonstrate this, we compare Gmc, 
the Monte-Garlo completeness for the MMF3 sample, using 
the cosmoOWLs profile variation prescription and effective 
beam variation, to Cerf, the semi-analytic ERF complete¬ 
ness. This comparison is shown in Fig. 7, where we show 
the difference between the two estimates as a function of 
Ysoo and ^500 as well the individual completeness values 
as functions of Ysoofor representative values of 05 oo slices 
through the 2D completeness and show the difference. 

The error function is a good approximation to the MG 
completeness for the cosmology sample, which uses a higher 
S/N cut and a larger Galactic mask than the full survey. 
The MC estimate corrects this analytic completeness by up 
to 20% for large resolved clusters, where Cmc is systemat¬ 
ically less complete than the ERF expectation, primarily 
due to variation in the cluster pressure profiles. For unre¬ 
solved clusters, the drop-off in Cmc is slightly wider than 
the ERF expectation, reflecting variation both of pressure 
profiles and of effective beams. 

The impact of these changes in completeness on ex¬ 
pected number counts and inferred cosmological param¬ 
eters for the cosmology sample is analysed in Planck 
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Collaboration XXIV (2015). The difference between the 
Monte-Carlo and ERF completeness results in a change in 
modelled number counts of typically ^ 2.5% (with a max¬ 
imum of 9%) in each redshift bin. This translates into a 
0.26(7 shift of the posterior peak for the implied linear fluc¬ 
tuation amplitude,(Tg. 

The MC completeness is systematically lower than the 
analytic approximation for the full survey. One of the causes 
of this is galactic dust contamination, which is stronger in 
the extra 20% of the sky included in the full survey area 
relative to the cosmology sample area. This tends to reduce 
the S/N of clusters on affected lines of sight. 

We note that this approach ignores other potential as- 
trophysical effects that could affect the completeness. Radio 
emission is known to be correlated with cluster positions, 
potentially ‘filling in’ the SZ decrement, though recent es¬ 
timates suggest that this effect is typically small in Planck 
data (Rodriguez-Gonzalvez et al. 2015). Departures of the 
pressure distribution from spherical symmetry may also af¬ 
fect the completeness, though this effect is only likely to 
be significat for nearby and dynamically disturbed clusters 
which may be large compared to the Planck beams. We 
test for some of these effects through external validation of 
the completeness in the next section, and explicitly through 
simulation in Sect. 4.5. 

Another source of bias is the presence of corre¬ 
lated IR emission from cluster member galaxies. Planck 
Collaboration XXIII (2015) show that IR point sources are 
more numerous in the direction of galaxy clusters, espe¬ 
cially at higher redshift, and contribute significantly to the 
cluster SED at the Planck frequencies. Initial tests, inject¬ 
ing clusters signal with the combined IR-|-tSZ spectrum of 
z > 0.22 clusters observed by Planck Collaboration XXIII 
(2015), suggest that this reduces the completeness for un¬ 
resolved clusters. Future work is warranted to characterise 
the evolution and scatter of this IR emission and to prop¬ 
agate the effect on completeness through to cosmological 
parameters. 

4.3. External validation of the completeness 

We validated our Monte-Carlo completeness calculation 
and our simple analytical ERF model by using the 
MCXC (Piffaretti et al. 2011) and SPT (Bleem et al. 2014) 
catalogues. The Planck detection threshold passes across 
the cluster distributions of these two samples. This is illus¬ 
trated in Fig. 4 of Chamballu et al. (2012) for the MCXC. 
This allows us to characterize our completeness by check¬ 
ing if the fraction of detected clusters follows the expected 
probability distribution as a function of their parameters. 
For each cluster of the MCXC catalogue, we use the MMF3 al¬ 
gorithm to extract its flux Ysoo and associated error cry at 
the location and for the size given in the x-ray catalogue. 
We then build the quantity (Ysoo — Wy)I t cty, q being 
the detection threshold (here 4.5) and cry the noise of the 
filtered maps. We make the corresponding histogram of this 
quantity for all the clusters and for the clusters detected by 
MMF3. The ratio of the two histograms is an empirical es¬ 
timate of the completeness. Results are shown In Fig. 8 
for the MCXC (left) and the SPT (right) catalogues. For 
MCXC, the estimation is in good agreement with the ex¬ 
pected simple analytical ERF model (0.5(1 -l-ERF)). For 
SPT, the estimated completeness is also in good agreement 
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Separation Carcminl 

Fig. 9. Cumulative distribution of angular separations be¬ 
tween estimated and input positions. The dashed vertical 
line denotes the Planck pixel size. 


except for (Ysoo — qcjY') jj cty > 1 where it is higher than 
the analytic expectation. We attribute this behaviour to the 
correlation between SPT and Planck detections. The SPT 
catalogue is a SZ-based, so a cluster detected by SPT will 
have a higher than random probability to be detected by 
Planck. This leads to an overestimation of the completeness 
at the high probability end. 

4.4. Position estimates 

We characterise the positional recovery of the Planck de¬ 
tections using injection into real data, including pressure 
profile and beam variation. We draw input clusters from a 
realistic distribution of (y 5 oo, 05 oo), the same as used for the 
reliability in Sect. 4.6. 

Fig. 9 shows the comparative performance of the indi¬ 
vidual detection codes, and of the reference position chosen 
for the union catalogue. PwS produces the most accurate 
positions, with 67% of detected positions being within 1.18 
arcmin of the input position. For MMFl and MMF3, the 67% 
bound is 1.58 arcmin and 1.52 arcmin respectively. The 
union and intersection accuracy follow that of the MMFs, 
with 67% bounds of 1.53 arcmin. We observe that our inter¬ 
code merging radius of 5 arcmin is conservative given the 
expected position uncertainties. 

4.5. Impact of cluster morphology 

Clusters are known to possess asymmetric morphologies 
and a wide range of dynamical states, from irregular merg¬ 
ing clusters to regular relaxed clusters. While the complete¬ 
ness simulations have included some morphology variations 
through variation of the injected radial pressure profile, this 
ignores the effects of the sub-structures and asymmetries, 
which may induce detection biases for large clusters at low 
redshift resolved by the Planck beams, FWHMRi 7 arcmin. 

Neither of the external samples used In Sect. 4.3 to val¬ 
idate the completeness allow us to properly probe resolved, 
irregular clusters at low-redshift. The MCXC is biased to¬ 
wards regular clusters due to X-ray selection effects and the 
Planck completeness drop-off lies substantially beneath the 
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Fig. 6. Completeness of the union and intersection samples at progressively lower S/N thresholds. From left to right, 
the thresholds are 8.5, 6.0 and 4.5 (the survey threshold). In the top panels, the dotted lines denote 15% completeness, 
the dashed lines 50% and the solid lines 85% completeness. In the bottom panels, the union is denoted by the diamonds 
with Monte-Carlo uncertainties based on binomial statistics, and the intersection is denoted by the solid lines. 


SPT mass limit at low redshift, so the drop-off is not sam¬ 
pled. 

We address the effects of realistic morphology by inject¬ 
ing into the Planck maps the raw 2D projected Compton-y 
signal from a sample of hydro-dynamically simulated cos- 
moOWLs clusters. The clusters were injected with a large 
enough angular extent, 6*500 = 20 arcmin, that they were re¬ 
solvable in the Planck data, and with a range of Ysoo that 
encompassed the expected completeness drop-off. 20 can¬ 
didate clusters were chosen from the sub-sample of cos- 
moOWLS clusters selected by the mass cuts discussed in 
Sect. 4.1 based on their dynamical state. The ten clusters in 
the sub-sample with highest kinetic-to-thermal energy ra¬ 
tio within 0500 constituted our disturbed sample, while the 
regular sample comprised the ten clusters with the lowest 
kinetic-to-thermal energy ratio within 05 oo- These clusters 
were injected 200 times, randomly distributed across the 
sky. We also created simulations injecting symmetric clus¬ 
ters with the UPP with the same parameters and locations 
as the hydro-dynamic projections. In all cases, the signals 
were convolved with Gaussian beams to separate the effects 
of beam asymmetries. 

The completeness for regular, disturbed and UPP clus¬ 
ters is shown for the union catalogue in the top panel of 
Fig. 10. There are no significant differences between the 
completeness functions for the regular and disturbed clus¬ 
ters. Both sets of hydro-dynamic clusters show a slight 
widening effect in the completeness caused by the variation 
in the effective pressure profile away from the UPP assumed 
for extraction (the same effect as discussed in Sect. 4.2). 

Morphology has a clear impact on the estimation of clus¬ 
ter position. The bottom panel of Fig. 10 shows the cumula¬ 
tive distribution of angular separation between union and 
input positions for the regular, disturbed and UPP clus¬ 
ters. The disturbed clusters show a significant reduction 


in positional accuracy. Part of this is physical in origin. 
The clusters centres were defined here as the position of 
the ’most-bound particle’, which traces the minimum of the 
gravitational potential and is almost always coincident with 
the brightest central galaxy. For merging clusters this posi¬ 
tion can be significantly offset from the centre of the peak of 
the SZ distribution. A matching radius of 10 arcmin, which 
is used in Sect. 7, ensured correct identification of detected 
and injected positions. 

4.6. Reliability 

The statistical reliability is the probability that a detection 
with given detection characteristics is a real cluster. We de¬ 
termine the reliability using simulations of the Planck data. 
Clusters are injected following the prescription in Sect. 4.1, 
except that the the clusters are injected such that cluster 
masses and redshifts are drawn from a Tinker et al. (2008) 
mass function and converted into the observable parame¬ 
ters (T 5 oo, 05 oo) using the Planck ESZ scaling 

relation (Planck Collaboration X 2011). The other com¬ 
ponents of the simulations are taken from FFP 8 simula¬ 
tion ensemble (Planck Collaboration XII 2015). The com¬ 
ponents include a model of galactic diffuse emission, with 
thermal dust (including some emission from cold-clumps), 
spinning dust, synchrotron and CO emission, and extra- 
galactic emission from the far infra-red background. The 
diffuse components are co-added to a set of Monte-Carlo 
realisations of CMB and instrumental noise. In addition to 
the cluster signal, we also inject point sources drawn from a 
multi-frequency model from the Planck sky model (Planck 
Collaboration XII 2015). These point sources are mock de¬ 
tected, using completeness information from the PCCS2 
(Planck Collaboration XXVI 2015), and harmonically in¬ 
filled using the same process as for the real maps prior to 
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Fig. 7. Differences between the semi-analytic and Monte-Carlo completenesses for MMF3. The left panels show the dif¬ 
ference for the full survey over 85% of the sky with a q = 4.5 threshold. The right panels show the difference for the 
MMF3 cosmology sample, covering 65% of the sky to a threshold of g = 6.0. The top panels show the difference MC — ERF 
in percent. The bottom panels compare the completenesses at particular 05 oo- the Monte-Carlo completeness is denoted 
by diamonds and the ERF completeness by solid lines. 




(^500,X-roy '^^Y(®5OO.X-roy))/SQRl'(2)/(7Y(05oo (t^soo.spl <1 *^y(®500,spi) )/SQRt(2)/c7Y(®500,spl) 

Fig. 8. MMF3 completeness for the PSZ2 catalogue (S/N threshold q > 4.5) determined from the MCXC (left) and SPT 
(right) catalogues. This external estimate (red histogram) is in good agreement with the analytic ERF calculation (solid 
blue line), except for SPT at the high probability end (see text). 
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Fig. 10. Impact of cluster morphology on the complete¬ 
ness (top panel) and position estimates (bottom panel) for 
resolved clusters. The simulated clusters are all injected 
with 0500 = 20', and all curves are for the union catalogue. 
Cluster morphology has no impact on the completeness, 
but a significant impact on the position estimation. 


SZ detection. This leaves a realistic population of resid¬ 
ual sources in the maps. After detection, candidates that 
lie within the simulated expanded source mask, or which 
match with the cold-core or IR source catalogues from the 
real data, have their S/N set to zero. 

Fig. 11 shows the reliability as a function of S/N for 
the union and intersection samples across the whole survey 
area and outside the 65% galactic mask used to define the 
cosmology samples. Relative to the PSZl, the reliability 
of the union has improved by 5%, the lower noise levels 
have revealed more real simulated clusters than spurious 
detections. As was the case with the PSZl, the reliability 
is improved significantly by removing more of the galactic 
plane, where diffuse and compact galactic emission induce 
extra spurious detections. 

4. 7. Neural network quality assessment 

We supplement the simulation-based reliability assessment 
with an a-posteriori assessment using an artificial neural- 
network. The construction, training and validation of the 
neural network is discussed fully in Aghanim et al. (2014). 



Fig.11. Cumulative reliability as a function of S/N. 


The network was trained on nominal mission Planck maps, 
with a training set composed of three elements: the po¬ 
sitions of confirmed clusters in the PSZl as examples of 
good cluster signal; the positions of PCCS IR and radio 
sources as examples of point-source induced detections; 
and random positions on the sky as examples of noise- 
induced detections. We provide for each detection a neu¬ 
ral network quality flag, Q_NEURAL = 1 — Qbad, following 
the definitions in Aghanim et al. (2014), who also tested 
the network on the unconfirmed detections in the PSZl. 
They showed that this flag definition separates the high 
quality detections from the low quality detections, as val¬ 
idated by the the PSZl external validation process, such 
that Q.NEURAL <0.4 identifies low-reliability detections 
with a high degree of success. 

459 of the 1961 raw detections possess QJ4EURAL < 
0.4 and may be considered low-reliability. This sample is 
highly correlated with the 1R_FLAG, with 294 detections in 
common. After removal of IR spurious candidates identified 
by the 1R_FLAG, as discussed in Sect. 3.2, we retain 171 
detections with bad Q_NEURAL, of which 28 are confirmed 
clusters. This leaves 143 unconfirmed detections considered 
likely to be spurious by the neural network. 

The QJ4EURAL flag is sensitive to IR induced spuri¬ 
ous: detections with low Q .NEURAL quality flag are clus¬ 
tered at low galactic latitudes and at low to intermediate 
S/N. This clustering is not seen for realisation-unique spu¬ 
rious detections in the reliability simulations, which are 
identifiable as noise induced. The reliability simulations 
underestimate the IR spurious populations relative to the 
QJ4EURAL flag. Gonversely, the neural network flag by 
construction does not target noise-induced spurious detec¬ 
tions: Qbad is the parameter trained to indicate IR-induced 
spurious. The neural network flag also has some sensitivity 
to the noise realisation and amplitude in the data: the as¬ 
sessment is different to that applied to the nominal mission 
maps in Planck Gollaboration XXXVI 2015. 

To place a lower limit on the catalogue reliability, 
we combine the QJMEURAL information with the noise- 
induced spurious detections from the reliability simulations. 
For each reliability simulation realisation, we remove the 
simulated IR spurious detections, which can be identified 
either as induced by the FFP8 dust component, and thus 
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Fig. 12. Lower limits on the catalogue reliability, esti¬ 
mated by combining the reliability simulations with the 
Q.NEURAL information (see text). 


present in multiple realisations, or as induced by injected 
IR point sources. We replace these spurious counts with 
the unconfirmed low Q_NEURAL counts, smoothed so as 
to remove the steps due to small number statistics. 

The combined lower limit of the reliability is shown in 
Eig. 12. The lower limit tracks the simulation reliability 
well outside the 65% galactic dust mask. Eor the whole 
survey region, the lower limit is typically 6 % lower than 
the simulation estimate, due either to over-sensitivity of 
the neural network to dusty foregrounds, or shortcomings 
in the FFP 8 galactic dust component. 

5. Parameter Estimates 

The SZ survey observable is the integrated Comptonisation 
parameter, Ygz- As was the case for the PSZl, each of the 
extraction codes has an associated parameter estimation 
code that evaluates, for each detection, the two dimensional 
posterior for the integrated Comptonisation within the ra¬ 
dius 5 i? 5 oo, L 5 R 500 , and the scale radius of the GNFW pres¬ 
sure 9s. The radius 5R500 is chosen as it provides nearly 
unbiased (to within a few percent) estimates of the total 
integrated Comptonisation, while being small enough that 
confusion effects from nearby structures are negligible. 

We provide these posteriors for each object and for each 
code, and also provide Torsoo in the union catalogue, de¬ 
fined as the expected value of the T 5 R 500 marginal distri¬ 
bution for the reference detection (the posterior from the 
code that supplied the union position and S/N). 

Below we also discuss the intricacies of converting the 
posteriors to the widely used X-ray parameters Tsoo—dsoo- 

5 . 1 . Y5R500 estimates 

To validate the I 5 R 500 estimates, we apply the posterior val¬ 
idation process introduced in Harrison et al. (2014) to the 
Tsrsoo marginal distributions. In brief, this process involves 
simulating clusters embedded in the Planck maps and eval¬ 
uating the Y — 9 posteriors for each (detected) injected clus¬ 
ter. For each posterior, we determine the posterior proba¬ 
bility, bounded by the contour on which the real under- 





30.0 

20.0 

15.0 

10.0 

8.0 

6.0 

4.5 


Fig. 13. Top panel: Results of the posterior validation 
for Ysrsoo- The histogram of the posterior probability, C,, 
bounded by the true Esrsoq parameter is almost uniformly 
distributed, except for a small excess in the tails of the 
posteriors, at C > 0.95. The histogram has been nor¬ 
malised by the expected counts in each C bin. Bottom panel: 
Comparison of recovered peak Ysrsoo to the injected Tsrsoo- 
The estimates are unbiased, though asymmetrically scat¬ 
tered, with a scatter that decreases as S/N increases. 


lying cluster parameters lie. If the posteriors are unbiased, 
the distribution of this bounded probability should be uni¬ 
formly distributed between zero and one. 

This process allows us to include several effects that vi¬ 
olate the assumptions of the statistical model used to esti¬ 
mate the posteriors. Firstly, by injecting into real sky maps, 
we include the non-Gaussian contributions to the noise 
on the multifrequency-matched-filtered maps that come 
from galactic diffuse foregrounds and residual point sources. 
Secondly, we include violations of the ‘signal’ model that 
come from discrepancies between the cluster pressure pro¬ 
file and the UPP assumed for parameter estimation, and 
from sky-varying and asymmetric effective beams that vary 
from the constant Gaussian beams assumed for estimation. 
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The clusters are injected using the process discussed in Sect. 

4.1, drawing injected pressure profiles from the set of cos- 
moOWLs simulated profiles. 

The top panel of Fig. 13 shows the histogram of C for the 
PwS Ysrsoo marginals. The distribution is flat, except for a 
small excess in the 0.95 — 1.0 bin, which indicates a small 
excess of outliers beyond the 95% confidence region, in this 
case 52% more than statistically expected. This suggests 
the posteriors are nearly unbiased, despite the real-world 
complications added to the simulations. Note that we have 
considered only posteriors where the injected FaRsoo > 
0.001 arcmin^, a cut that removes the population effects of 
Eddington bias from consideration: we focus here on the 
robustness of the underlying cluster model. 

The bottom panel of Fig. 13 shows the peak recovery 
from the PwS F 5 R 500 marginals compared to the true in¬ 
jected values. The peak estimates are unbiased relative to 
the injected parameters. 


5.2. Conversion to Ysoo 

The (F 5 R 5 oo, 0 s) estimates can be converted into (Fsooj^soo) 
estimates using conversion coefficients derived from the 
UPP model that was assumed for extraction. However, 
when the underlying pressure distribution deviates from 
this model, the conversion is no longer guaranteed to accu¬ 
rately recover the underlying (Fsooj^soo) parameters: vari¬ 
ation of the pressure profile can induce extra scatter and 
bias in the extrapolation. 

We demonstrate this by applying the posterior val¬ 
idation process to the ¥500 posteriors, defined as the 
^ 5 R 500 posteriors scaled with the UPP conversion coeffi¬ 
cient, as estimated from injected clusters whose pressure 
profiles are drawn from the cosmoOWLS pressure profile 
ensemble. We validate posteriors for ¥500 calculated in two 
ways: firstly by marginalising over the 6500 parameter, re¬ 
ferred to in previous publications as ‘Y blind’; and secondly 
by slicing the (¥ 500 ,^ 500 ) posteriors at the true value of 
^500 j equivalent to applying an accurate, externally mea¬ 
sured delta-function radius prior. 

Fig. 14 shows the bounded probability histograms for 
the two Y 500 posteriors and Fig. 15 shows the scatter of the 
peak of the posteriors with the input values of ¥500- The 
marginal Y 500 posteriors are poor, with histograms skewed 
towards the tails of the distribution and large numbers of 
> 2(7 outliers. The scatter plot reveals the peak estimates to 
possess a large scatter and to be systematically biased high. 
In contrast, the peak p(y 5 oo| 05 oo) estimates have much bet¬ 
ter accuracy and precision and are distributed around the 
input values with low scatter. The bounded probability his¬ 
togram ofp(¥ 5 oo| 6 * 5 oo) shows that while there is a noticeable 
excess of detections in the wings, the posteriors are reason¬ 
ably robust. If the posteriors were Gaussian, the skewness 
of the p(Y 5 oo 1 ^ 500 ) histogram towards the tails would be 
consistent with an underestimate of the Gaussian standard 
deviation of 21 % 

We therefore recommend that, to estimate ¥500 accu¬ 
rately from Planck posteriors, prior information be used 
to break the (¥500,6*500) degeneracy. However, we note that 
the uncertainties on such Y 500 estimates will be slightly un¬ 
derestimated. 





Fig. 14. Bounded probability histograms, as in the top 
panel of Fig. 13, but for the converted p(y5oo) marginal 
and p(¥50016*500) sliced posteriors. 


5.3. Mass and ¥500 estimates using scaling priors 

The key quantity which can be derived from SZ observables 
is the total mass of the detected clusters within a given 
overdensity (we used A = 500). To calculate the mass from 
Planck data it is necessary to break the size-flux degeneracy 
by providing prior information, as outlined in the previous 
section. We used an approach based on Arnaud et al. (in 
prep.), where the prior information is an expected function 
relating ¥500 to @500 that we intersect with the posterior 
contours. We obtained this relation by combining the defi¬ 
nition of M 500 (see Eq. 9 in Planck Collaboration XX 2014, 
connecting M 500 to 6 * 500 , for a given redshift z) with the 
scaling relation ¥500 — M 500 found in Planck Collaboration 
XX (2014). A similar approach was also used in Planck 
Collaboration XXIX (2014), but in this work we use the full 
posterior contours to associate errors to the mass value. 

We illustrate our method in Fig. 16. At any fixed value 
of 6 * 5 , we study the probability distribution and derive the 
dsRsoo associated to the maximum probability, i.e. the ridge 
line of the contours (red continous line in Fig. 16). We also 
derive the ¥ 5^500 limits enclosing a 68 % probability and 
use them to define a upper and lower degeneracy curve 
(dashed lines). From the intersection of these three curves 
with the expected function (cyan line), we derive the Msz 
estimate and its Icr errors, by converting Y 5 _r 5 oo to F 500 and 
then applying the F 500 — scaling relation prior at the 
redshift of the counterpart. 

Msz can be viewed as the hydrostatic mass expected 
for a cluster consistent with the assumed scaling relation, 
at a given redshift and given the Planck {Y — 9) posterior 
information. We find that this measure agrees with exter¬ 
nal X-ray and optical data with low scatter (see Sect. 7). 
For each Msz measurement, the corresponding Y 500 from 
the scaling relation prior can be calculated by applying the 
relation. 

We underline that the errors bars calculated from this 
method consider only the statistical uncertainties in the 
contours, not the uncertainties on the pressure profile nor 
the errors and scatter in the ¥500 — M scaling relation, 
and should thus be considered a lower limit to the real 
uncertainties on the mass. 
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Fig. 15. Scatter of the recovered estimates of Ysoo with 
the input Ysoo- Top panel: for the marginalised Ysoo pos¬ 
terior, ‘Y blind’. Bottom panel: for the sliced posterior 
p(^ 50 o|^ 50 o)) assuming an accurate radius prior. 


We used the masses for the confirmation of candidate 
counterparts (see Sect. 7) and we provide them, along with 
their errors, in the PSZ2 catalogue for all detections with 
confirmed redshift. We compared them with the masses pro¬ 
vided in PSZl for the detections where the associated coun¬ 
terpart (and thus the redshift value) has not changed in the 
new release (see Appendix B). We find very good agreement 
between the two values which are consistent within the er¬ 
ror bars over the whole mass range. 

In the individual catalogues, we provide for all entries 
an array of masses as a function of redshift (Msz(.z)), which 
we obtained by intersecting the degeneracy curves with the 
expected function for different redshift values, from z = 0 
to z = 1. The aim of this function is to provide a useful 
tool for counterpart searches: once a candidate counterpart 



Fig. 16. Illustration of the posterior probability contours 
in the Ys^soo ~ plane for a cluster detected by Planck: 
the contours show the 68, 95 and 99 percent confidence 
levels. The red continuous line shows the ridge line of the 
contours while the dashed lines the l-cr probability value at 
each Og. The cyan line is the expected Y — 9 relation at a 
given redshift that we use to break the degeneracy. 


Table 3. Results of fits between S/N from the PSZl and 
PSZ2, using the fitting function in Eq. 6. The assumed cor¬ 
relation of the uncertainties of si and S 2 was 0.72. 


Si 

S2 

A a a 

PSZ2 

PSZl 

0.76 ±0.08 0.72 ±0.01 0.53 ± 0.02 


is identified, it is sufficient to interpolate the Msz(z) curve 
at the counterpart redshift to estimate its mass. 

6. Consistency with the PSZl 

The extra data available in the construction of the PSZ2 im¬ 
proves the detection S/N and reduces statistical errors in 
the parameter and location estimates. Here we assess the 
consistency between the two catalogues, given the matching 
scheme discussed in Sect. 7.1. 


6.1. Signal-to-noise 

We fit the relation between S/N for common PSZl and 
PSZ2 using the the approach and model discussed in Sect. 
3.4. For the PSZl and PSZ2, the likelihoods for si and 
S 2 have a strong covariance, as more than half of the 
PSZ2 observations were used in the construction of the 
PSZl. We therefore assign a covariance of 0.72 between 
the two S/N estimates, as is appropriate for Gaussian er¬ 
rors sharing 53% of the data. As the errors are not truly 
Gaussian, we allow for an intrinsic scatter between the 
S/N estimates to encapsulate any un-modelled component 
of the S/N fluctuation. 

The consistency of the S/N estimates between the 
PSZl and PSZ2 are shown in Fig. 17 and the best fit¬ 
ting model in is shown in Table 3. Detections with PSZ2 
S/N > 20 are affected by changes in the MMF3 S/N defi¬ 
nition. For the PSZl, the empirical standard deviation of 
the filtered patches was used to define the S/N in this 
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Fig. 17. Comparison of S/N values for common PSZl and 
PSZ2 detections. The best fit relation is plotted in red, with 
2cr scatter plotted by dashed red lines. The green dashed 
line denotes the 1-1 relation. 


regime, while the theoretical standard deviation of Melin 
et al. (2006b) was used for lower S/N. MMF3 now uses the 
theoretical standard deviation for all S/N, consistent with 
the ESZ and the definitions in the other detection codes. 
For this reason, the best fit model ignores detections at 
S/N > 20 in either catalogue. The MMF3 S/N show a flat 
improvement relative to the ESZ S/N (which was produced 
solely by MMF3), consistent with the reduced noise in the 
maps. 

If the Compton-Y errors are entirely Gaussian in their 
behaviour, we should expect the S/N to increase by 37% 
between the PSZl and PSZ2, ie: a = 0.73. This is consistent 
within Icr with the fit, which describes the S/N behaviour 
well to S/N< 20. 

6.2. Position estimation 

The distribution of angular separations between the 
PSZ2 and PSZl position estimates is shown in Fig. 18. 
Of the common detections, 80% of the PSZ2 positions he 
within one Planck map pixel width, 1.7 arcmin, of the 
PSZl position. MMF3 does not allow for sub-pixel position¬ 
ing, so if the MMF3 position was used for the union in both 
the PSZl or PSZ2, the angular separation will be a mul¬ 
tiple of the pixel width. This is evident in the cumulative 
distribution of angular separations as discontinuities at 0, 
1 and y/2 pixel widths. 

We also compare the position discrepancy between 
the SZ detection and the X-ray centres from the MCXC 
(Piffaretti et al. 2011). The bottom panel of Fig. 18 shows 
the distributions of these angular separations for the 
PSZ2 and PSZl. The distributions are calculated from the 
full MCXC match for each catalogue: the PSZ2 includes 
124 new detections. The PSZ2 position estimates are clearly 
closer to the X-ray centres than the PSZl: for the PSZl, 
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Fig. 18. Top panel: Separation between PSZ2 and 
PSZl positions for common detections. Middle panel: 
Cumulative distribution of the angular separation between 
PSZl and PSZ2 positions, with the Planck pixel width in¬ 
dicated by a dashed vertical line. Bottom panel Cumulative 
distributions of angular separation to MCXC x-ray centres, 
for all PSZl and PSZ2 MCXC matches. The vertical dashed 
line denotes the Planck Healpix pixel size. 


the 67% error radius is 1.85 arcmin, while for the PSZ2 this 
reduces to 1.6 arcmin. 
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6.3. Missing PSZl detections 

The PSZl produced 1227 union detections. While the num¬ 
bers of detections has increased by 35% in the PSZ2 to 
1653, the number of common detections is 936: 291 (23.7%) 
of the PSZl detections disappear. The high-purity intersec¬ 
tion sample loses 44 detections, of which 20 are lost entirely, 
and 24 drop out of the intersection after one or two codes 
failed to detect them. In this section, we discuss these miss¬ 
ing detections. Table E.l details each of the missing detec¬ 
tions and provides an explanation for why each is missing. 

The first type of missing detection are those that fall un¬ 
der the new survey mask, due to the increase in the number 
of point sources being masked. The masked areas sare pre- 
processed with harmonic infilling to prevent spurious detec¬ 
tions induced by Fourier ringing. The increase of the mask 
area is driven by S/N improvements for IR sources in the 
high frequency channels. While the increase in the masked 
area is small (0.1% of the sky), the correlation between IR 
point sources and galaxy clusters leads to a larger percent¬ 
age of detections being masked. In the PSZl, these detec¬ 
tions were contaminated by point source emission, but the 
emission was just beneath the point source masking thresh¬ 
old. 21 PSZl union detections fall behind the new mask. Of 
these, three were confirmed clusters, none received the high¬ 
est validation quality flag of 1 (denoting probable clusters) 
in the PSZl validation process^, four received the interme¬ 
diate validation quality flag 2, and 14 received the lowest 
validation quality flag of 3, denoting probable spurious. 

The second type of missing detection is one which has a 
matching detection in the full-mission data, but where the 
detection was rejected either by the infra-red spurious cuts 
or by PwS internal consistency cuts, both of which are dis¬ 
cussed in Sect. 3.2. The IR cuts are responsible for cutting 
33 unconfirmed PSZl detections, of which six were in the 
intersection sample. In the PSZl validation process, none 
of these received validation quality flag of 1, seven received 
quality flag 2 and 26 received quality flag 3. These were 
all S/N < 7 detections. Five detections were lost because 
PwS was the only detecting code in the PSZ2 and they failed 
PwS consistency criteria: two of these were confirmed clus¬ 
ters. 

The final type contains the majority (232) of the miss¬ 
ing detections. These are low-significance detections close 
to the PSZl threshold that have downward-fluctuated with 
the full mission data and are now beneath the PSZ2 thresh¬ 
old. This occurs for some detections despite the fact that 
the S/N improves for most. The top panel of Fig. 19 shows 
the PSZl S/N distribution of the downward-fluctuated de¬ 
tections. These were weak detections: 87% were within 0.5 (t 
of the detection threshold and 82% of them were single¬ 
code detections. While many of these may be spurious de¬ 
tections, 81 confirmed clusters have been lost. 61 of these 
were single-code detections and 70 of them were within 0.5 (t 
of the threshold. Based on Planck data alone, these clusters 
were weak SZ detections and were likely to be Fddington 
biased above the threshold in the PSZl. We have estimated 


^ The PSZl validation process produced three quality flags for 
unconfirmed clusters. These were based on a combination of SZ 
signal quality, X-ray signal in the RASS maps and IR signal in 
the WISE maps. Class 1 candidates satisfied good quality in all 
three measures and were high reliability candidates. Class 2 sat¬ 
isfied at least one measure with good quality, while class 3 failed 
all three measures and so were considered probably spurious. 
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Fig. 19. Distribution of S/N for missing nominal mission 
detections lost due to downward fluctuation of the S/N 
rather than because of spurious rejection cuts or changes in 
the survey mask. Top panel Detections lost from the PSZl. 
Bottom panel Detections lost in simulations of the transi¬ 
tion from the nominal to full mission. 


the S/N for these lost PSZl detections in the full-mission 
maps using PwS in a non-blind analysis at the PSZl po¬ 
sitions. Fig. 20 shows the distribution of these non-blind 
S/N: for most an apparently significant signal still exists 
in the maps, but it is now too weak to exceed the detec¬ 
tion threshold, typically lying between 2 < S/N < 4. The 
non-blind S/N for this category is shown per detection in 
Table E.l. Two detections have a non-blind S/N above the 
selection threshold. For these detections, the noise level for 
the non-blind analysis (centred on the PSZl location) was 
lower than for any of the patches in the mosaic used for the 
cluster detection. 

To verify that this sample of missing detections is con¬ 
sistent with the change in data, we simulated the transi¬ 
tion from the PSZl to the PSZ2 using FFP8 half-mission 
noise realisations to approximate the nominal mission: this 
produces a pair of data-sets with appropriately correlated 
noise characteristics. A common sample of clusters and 
point sources were injected into the simulations and the 
full pipeline was applied to construct catalogues from both 
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S/N non-blind 

Fig. 20. Non-blind PwS S/N for the SZ signal at the location 
of missing PSZl detections that were not masked out or cut 
for IR contamination. 


simulated datasets. The simulations produced a total loss 
of 353 detections, of which 24 were lost due to the expan¬ 
sion of the point source mask, ten were lost due to changes 
to the PwS spurious rejection criteria, and 319 were lost due 
to downward fluctuation of the S/N beneath the detection 
threshold. 

The S/N distribution of the latter group is shown in the 
bottom panel of Fig. 19. 75% lie within 0.5 ct of the detection 
threshold and 85% were single-code detections. While this 
group was primarily composed of 230 spurious detections, 
89 injected clusters were lost. The loss of these injected 
clusters illustrates that, as a statistical process, cluster de¬ 
tection is dependent on the realisation of the noise in the 
filtered patch-maps and we should expect that substantial 
numbers of confirmed but weak cluster detections will be 
lost due to noise fluctuations. 

These simulations over-estimate the loss-rate of nominal 
mission detections. This may be due in part to unsimulated 
changes in the sample selection applied to the real data. 
We were unable to simulate systematic changes in the IR 
spurious rejection that may, had they been incorporated, 
have resulted in some spurious detections from the nominal 
mission simulation being cut from the comparison. 

6.4. Compton Y estimates 

The Compton Ysrsoo estimates from each code are com¬ 
pared to the PSZl estimates in Fig. 21. The Fsrsoo esti¬ 
mates that we consider here are the mean estimates of the 
ksRsoo marginal posteriors, having marginalised over the 
scale radius Og. 

The best fit relations between the PSZ2 and PSZl values 
for each code are shown in Table 4. These were fit using a 
similar procedure to the S/N estimates discussed in the 
previous section. We assume a log-linear relation between 
the estimates of the form 

\og^ = A + a\og^, (7) 


with a log-normal intrinsic scatter dint and Tpiv = 3 x 10 ^ 
arcmin^. We again assume a bivariate Gaussian likelihood 
for the estimates, with a correlation of 0.72. 

Fig. 21 compares the Tsrsoo estimates for each of the 
three detection codes. High S/N detections have more con¬ 
sistent estimates of Ysrsoo- For MMF3, detections at S/N> 20 
are significantly changed due to the changes in the treat¬ 
ment of these detections discussed in Appendix C. These 
points are excluded from the fit to the relation. The scatter 
on the high S/N estimates is determined by the robustness 
of the noise power spectrum estimation to small changes in 
the data. For PwS, the high S/N estimates have particularly 
low scatter, due to the robust nature of the noise estima¬ 
tion that accounts for Compton-Y ‘noise’ contributed by 
neighbouring clusters. 

The low S/N detections show systematic deviations for 
each of the codes. For the MMFs, these are caused by the 
correction of PSZl Eddington bias in the PSZ2 data, which 
is visible in Fig. 21 as clouds of faint points where the 
ksRsoo estimate reduces in the PSZ2. The opposite is the 
case for PwS estimates, where the faint detections show up¬ 
ward deviation in the PSZ2. This is caused by a change 
in the priors: for the PSZl, PwS used a power-law prior in 
Y 5 R 500 , which was replaced in the PSZ2 with the uninfor¬ 
mative flat prior, as this produced more robust T 5 R 500 es¬ 
timates in the posterior validation process discussed in the 
previous section. We have confirmed that PwS behaves in 
the same way as the MMFs when uninformative priors are 
used for both PSZl and PSZ2 parameter estimates (see the 
bottom right panel of Fig. 21). 

To verify that the bias effects seen in Fig. 21 are within 
expectations, we extracted Ysrsoo estimates from the half 
to full-mission transition simulations described in Sec. 6.3. 
We confirm the same behaviour in these simulations as in 
the real data: low S/N detections from the MMFs show 
a correction of Eddington bias in the full mission while 
PwS low S/N detections are affected by change from power- 
law to uninformative priors in the posterior estimation and 
are typically higher in the full mission. 


7. Ancillary Information 

7.1. Cross-match with PSZl 

We begin the search for counterparts by conducting a cross¬ 
match with the well-validated PSZl. All matches within 5 
arcmin of a PSZl detection are accepted as a true match. 
Both catalogues used this radius as the merging limit to 
define unique detections, both in the merge of Cartesian 
patch catalogues to form an all-sky catalogue and in the 
formation of the union). This step produced no multiple 
matches. 

Several of our detections are clear matches with 
PSZl detections at higher radii than this, so we consider 
matches out to 10 arcmin, as is the case with the X-ray 
and optical counterpart searches described below. This step 
produced 18 potential matches, two of which were non¬ 
unique. We apply a further condition to accept these high- 
separation matches: that the PSZ2 S/N be greater than the 
PSZl S/N and be consistent with the S/N relation deter¬ 
mined in Sect. 6.1. For the two non-unique matches, the 
nearer match was chosen both times and this match also 
better fit the S/N relation. 
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Fig. 21. Comparison of I 5 R 500 estimates from individual codes in the PSZl and PSZ2. The I 5 R 500 estimator is the mean 
of the Y — 9 posteriors, marginalised over 6s (‘Y blind’). The circled red points denote sources with S/N > 20. The 
dashed green lines show the 1 — cr envelope of the best-fit relations shown in Table 4. MMFl estimates are shown top left, 
MMF3 top right. PwS estimates are shown bottom left. The bottom right panel compares PwS estimates having re-analysed 
PSZl data using uninformative priors on Ysrsoo and 0s- 


Table 4. Results of fits between Y 5 R 500 from the PSZl and PSZ2, following Eq. 7. 


Yi 

Y 2 

A 

a 

Cint 

PSZl MMFl 

PSZ2 MMFl 

-0.087 ±0.006 

1.00 ± 0.02 

0.083 ± 0.003 

PSZl MMF3 

PSZ2 MMF3 

-0.054 ±0.002 

1.05 ± 0.01 

0.054 ± 0.003 

PSZlPwS 

PSZ2 PwS 

0.056 ±0.005 

1.02 ±0.01 

0.068 ± 0.003 


7.2. X-ray information 

We use the Meta-Catalogue of X-ray deteeted Clusters of 
galaxies (MCXC, Piffaretti et al. 2011) for the association 
of Planck SZ candidates with known X-ray clusters, as 
was done in Planck Collaboration XXIX (2014). MCXC is 
based on the ROSAT All Sky Survey and complemented 
with other serendipitous catalogues and with the Einstein 
Medium Sensitivity Survey. It includes 1743 clusters 
distributed over the whole sky and provides coordinates, 
redshifts and X-ray luminosity measured within R 500 , 
-^x, 500 - The association of Planck SZ candidates with 
MCXC clusters follows two steps: first a positional match¬ 
ing between the catalogues then a verification of the 
association using the Tx.soo ~ relation (Pratt et al. 


2009). In the first step, we looked for possible counterparts 
of Planck SZ candidates in the MCXC within a searching 
radius of 10 arcmin around the Planck position. We 
found one counterpart for 537 candidates and multiple 
matches for another 16 objects. In the second step, we 
verified our associations by looking at their position in the 
Lx ,500 — Af^oo plane (Fig. 22). For the X-ray luminosity, 
we use the Ex ,500 value provided in the MCXC, while 
we calculate the mass from our own data, as described 
in Sect. 5. In Fig. 22, we compare our results with the 
expected Tx.soo — Yfsoo relation (Piffaretti et al. 2011): 
we consider as good associations those whose position in 
the Lx ,500 — does not differ from the expected one 

by more than twice the intrinsic scatter in the relation 
(^int = 0.183 Pratt et al. 2009). Based on this criterion we 
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discarded the association with an MCXC cluster for two 
objects, PSZ2 G086.28 + 74.76 and PSZ2 G355.22 - 70.03, 
both new PSZ2 detections. 

We performed a further check of the candidate coun¬ 
terparts, by studying the separation between the Planck 
and the MGXG positions. Indeed, the relatively large 
search radius (10 arcmin) may have led to spurious 
associations, which might have escaped our selection on 
the Lx ,500 — relation. In Fig. 23, we compare the 

separation between the Planck and the MGXG positions 
with two relevant angular scales: the positional uncertainty 
of the Planck detections 0err (90% confidence level, pro¬ 
vided in the catalogue) and the cluster size as quantified 
by 6 * 500 ®- Ideally, one would keep as good counterparts 
those systems where the angular separation is smaller than 
both 0500 and 0err (lower left quadrant in Fig. 23) but this 
would lead to a large number of rejected matches including 
many objects in the PSZl. Therefore, we chose a less 
conservative criterion: we excluded only those associations 
where the separation is larger than both 05 oo and 0 err 
(upper right quadrant in Fig. 23). We thus allow the 
separation to exceed 05 oo if the MGXG counterpart falls 
within the Planck accuracy (upper left quadrant) and to 
exceed Oerr if it is smaller than the cluster expected size 
(lower right quadrant). We noticed that the most deviant 
clusters in the latter case (with 6 > 205 oo) are associated 
with nearby clusters {z < 0.14) with 05 oo > 7 arcmin, 
resolved by Planck. In this phase of the analysis, we thus 
discarded the three associations to PSZ2 G247.97 -I- 33.52, 
PSZ2 G212.93 - 54.04, and PSZ2 G209.79 -b 10.23 in the 
upper right quadrant of Fig. 23: for these systems, the 
separations between the Planck and MGXG are always 
larger than 5 arcmin, which would correspond to a physical 
distance of ~ IMpc at the redshift of the MGXG objects 
(all at z > 0 . 2 ). 

We also used the position in the Tx.soo ~ -^soo and in 
the separation plane (Fig. 23) to select the most likely coun¬ 
terparts for the objects where two or more MGXG clusters 
were found within our search radius of 10 arcmin. For seven 
out of 16 objects, one counterpart does not match the cri¬ 
teria described above and we are thus left with only one 
good counterpart. For six Planck detections, both MCXC 
counterparts fullfill our requirements. We thus rank them 
based on their distance from the Tx ,500 ~ -^500 scaling re¬ 
lation and their separation in terms of 05 oo and 0 err, and 
select as most likely counterpart the one with smaller values 
for at least two out of three indicators. We provide details 
of the other possible counterparts in the Comments file. In 
two cases, the same MCXC cluster can be associated with 
two Planck detections and we used the procedure described 
above to select the most likely associations. 

In the last step of our analysis, we checked our matching 
with MCXC with the matches made in the PSZl catalogue: 
in most cases the MCXC counterparts in the two catalogues 
coincide. We examined in detail the cases where, following 
our selection criteria, we would have broken the association 
with the MCXC counterpart which was chosen in PSZl. 
This led to three restorations: 


® We calculated 05 oo from the mass proxy M 500 , using the 
redshift of the MCXC counterpart. 


- PSZ2 G247.97-b33.52 (PSZl index 842) association with 
RXCJ0956.4 — 1004 lies in the forbidden area in the 
separation plane. However, RXCJ0956.4 — 1004 (also 
known as A901) is a multi-component cluster (Bosch 
et al. 2013), and the PSZ2 position lies close to the po¬ 
sition of one of the components. We thus decided to keep 
the association. 

- PSZ2 G302.41 -b 21.60(PSZ1 index 1054) and PSZ2 
G332.29— 23.57(PSZ1 index 1158) are both associated 
very low-redshift clusters (RXC J1248.7—4118 and RXC 
J1847.3 —6320, both z < 0.015) which are marginal out¬ 
liers in the L-M plane. However, our mass proxy esti¬ 
mate may be less reliable for local objects due to the 
large cluster extent and we thus decided to keep these 
associations. 


7.2.1. Comparison to L — M relation 

It is interesting to note in Fig. 22 that most points lie be¬ 
low the expected scaling relations of Pratt et al. (2009), al¬ 
though well within the intrinsic scatter, meaning that clus¬ 
ters in our subsample are systematically under-luminous 
(by about 21%, or —0.41tTint) in X-rays at a given mass. We 
recall here that this subsample (intersection of PSZ2 with 
MCXC) is not representative and thus cannot be quantita¬ 
tively compared with a well defined representative sample 
such as REXCESS, for which the Tx.soo — -^soo was de¬ 
rived by Pratt et al. (2009). The systematic offset observed 
in Fig. 22 does not contradict the good agreement between 
X-ray predictions and Planck measurements found with a 
statistical approach in Planck Collaboration X (2011). It 
can be explained taking into account selection effects and 
the scatter in the Y — Lx scaling relation: when cutting in 
SZ S/N, clusters with a high SZ signal (and thus a high 
mass) for a given X-ray luminosity are preferentially se¬ 
lected. Another effect, which could also partly contribute 
to the offset in Fig. 22, is the presence of a cluster popu¬ 
lation with different X-ray properties in the PSZ2 sample 
which will be discussed in Sect. 8.2. 


7.3. Optical information 

We benefit from a wealth of publicly available data over 
the northern sky, principally thanks to the Sloan Digital 
Sky Survey (SDSS, York et al. 2000) that covers most of 
the northern extragalactic sky with imaging in five optical 
bands (ugriz). A number of cluster catalogues have been ex¬ 
tracted from these data using different finding algorithms 
(Koester et al. 2007; Hao et al. 2010; Wen et al. 2012). 
Among these, the redMaPPer catalogue (Rykoff et al. 2014), 
published since the PSZl release and containing many more 
clusters, has proven to be the most useful for identifying 
counterparts to Planck SZ sources. We also supplement 
redMaPPer with other optically information. 

7.3.1. RedMAPPer 

The redMaPPer algorithm detects clusters by looking for 
spatial over-densities of red-sequence galaxies. It provides 
accurate photometric redshift estimates for all sources, 
spectroscopic redshifts for the brightest central galaxy 
(BCG) when available, and richness estimates. We used 
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Fig. 22. Comparison of candidates associated with the 
MCXC catalogue with the expected Lx ,500 — -^500 scaling 
relation (red line). The parallel dashed lines identify the 
region of the plane within 2(Tint from the expected scaling 
relation, where aint is the logarithmic intrinsic scatter of 
the relation we used. Black points are confirmed MCXC 
associations, while magenta squares mark the associations 
discarded by the L-M criterion. Pairs of coloured diamonds 
mark the two possible counterparts for objects with multi¬ 
ple associations. 



Fig. 23. Separation between the Planck and MCXC posi¬ 
tions in terms of the positional uncertainty of the Planck 
detection and of the cluster i? 5 oo- The horizontal and ver¬ 
tical dashed lines mark our acceptance threshold. 


the proprietary redMaPPer catalogue (v5.10) provided by 
the authors and containing over 400,000 objects. 

In our procedure, detailed further in Bartlett et al. 2014 
(in prep.), each Planck SZ source is first matched to a max¬ 
imum of three redMaPPer clusters falling within a radius of 
10 arcmin. They are subsequently ranked by richness and 
labeled first-, second-, and third-ranked matches. We then 
calculate the Planck mass proxy, Mgz, for each SZ source at 
the redshifts of its matched redMaPPer clusters. The best 
redMaPPer counterpart is then selected based on cuts in an¬ 
gular separation and richness. The angular cuts incorporate 
both the Planck positional uncertainty and the physical ex¬ 
tent of the cluster estimated from the calculated M^z- 

These angular criteria alone would leave multiple pos¬ 
sible counterparts in many cases, given the high surface 



D InA/oin 

Fig. 24. Distribution of positional matches within a 10 ar¬ 
cmin radius in the richness-Mgz plane. The red points in the 
upper panel represent the highest richness match, blue the 
second (when present) and green points the third-ranked 
richness match (when present). The mean scaling law from 
Rozo et al. (2014) is shown as the solid line, with the dashed 
lines delineating the ±3cr band. In the lower panel, we show 
the distribution of these points relative to the mean rela¬ 
tion, normalized to the logarithmic scatter. The red, blue 
and green histograms refer to the first-, second-, and third- 
ranked matches, respectively. 


density of redMaPPer clusters. Any ambiguity is efficiently 
reduced by the richness cut, which is based on the existence 
of a well-defined relation between richness. A, and Msz- The 
relation was established by Rozo et al. (2014) on the Planck 
2013 SZ cluster catalogue and is expressed as 

{\nX\Msz) = a + (8) 

with a = 4.572 ± 0.021, a = 0.965 ± 0.067, and Mp = 
5.23 X IO^'^Mq. The measured dispersion at given Msz is 
CTin A = 0.266 ±0.017. 

In Fig. 24 we compare the distribution of the first, 
second-, and third-ranked redMaPPer matches in the A — 
Msz plane to this scaling relation. The quantity Ain A = 
[ln(A) — (lnA)]/crinA is the deviation of measured rich¬ 
ness from the expected mean. We see that the hrst-ranked 
matches (in red) display a prominent peak with ±3(7inA- 
This reaffirms the existence of the scaling relation and mo- 
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Fig. 25. Selection criteria for redMaPPer counterparts. 
Objects are plotted following the same colour scheme as 
in the previous figure, i.e., first-, second- and third-ranked 
matches represented by red, blue and green symbols, re¬ 
spectively. The bands delineate the acceptable region de¬ 
fined by the angular criterion, and circled points indicate 
objects that also satisfy the richness cut (Criterion 1). 
Upper panel - Single matches within a 10 arcmin radius. 
Lower panel - Double matches within a 10 arcmin radius. 
Note that in both panels, only first-ranked matches appear 
as good counterparts by satisfying both criteria. 


tivates its use in defining the final redMaPPer counterparts 
for the PSZ2. 

We define the best redMaPPer counterparts with the 
following cuts: 

1. (0/0err) < 1 OR (0/05Oo) < 1; 

2. |AinA| < 3; 

3. When more than one object remains a possible counter¬ 
part, we choose the highest ranked match. 

The criteria on angular separation, 9, allow objects with 
centres either within Planck's, positional uncertainty or 
within the estimated size of the cluster (or both). The sec¬ 
ond criterion imposes the richness requirement based on 
the scaling relation. If there remain more than one possible 
counterpart satisfying these two criteria, then we choose 
the one with the largest richness. These latter cases, how¬ 
ever, deserve closer examination, in particular for potential 
projection effects. 


Table 5. Distribution of the 375 good counterparts. 


Matches 


Single Double Triple 


Rank (/40) (/85) (/438) 


1 . 17 58 0 283 8 0 

2 . 0 6 0 

3. 2 


One object with three possible good matches. 


Fig. 25 shows the distribution of redMaPPer objects for 
Planck sources with single and double matches. The dashed 
lines delineate the angular cuts, while open circles identify 
those objects that also satisfy the richness cut. We see that 
the richness cut effectively eliminates objects that would be 
accepted on angular criteria alone. For the double matches, 
there are no good second-ranked redMaPPer counterparts 
because only the highest richness object satisfies the rich¬ 
ness cut. 

Table 5 summarizes the distribution of the 375 counter¬ 
parts found with the above cuts. It is grouped into sets of 
columns for single-, double-, and triple-matched Planck ob¬ 
jects. For each grouping, the total number of Planck sources 
is given in parentheses in the heading. The table entries list 
the number of matches in matrix format as follows: the el¬ 
ement (i, j) of a matrix gives the number of Planck sources 
with both i- and j-ranked redMaPPer counterparts; for ex¬ 
ample, six of the 438 triple matches have a good second- 
ranked counterpart only, while eight have both good first- 
and second-ranked counterparts. 

7.3.2. Other optical information 

We perform targeted searches for counterparts within 
the SDSS footprint for all Planck sources without good 
redMaPPer matches by applying the redMaPPer algorithm 
on a case-by-case basis. This yielded an additional 17 coun¬ 
terparts and associated redshifts. 

We add optical confirmations and redshifts from sev¬ 
eral other sources. These include optical counterparts for 
PSZl clusters published recently from PanSTARRs (Liu 
et al. 2014) and from Planck collaboration optical follow¬ 
up observations (Planck Collaboration et al. 2014; Planck 
Collaboration Int. XXXV 2015). We also search for coun¬ 
terparts in the NED® database, again removing any negated 
duplicate matches from the systematic searches. We com¬ 
pared the NED redshifts for PSZ2 matches to the redshifts 
from all the other ancillary catalogues we have studied. 
The NED redshifts have 88% agreement with these sources 
within Az < 0.02. We therefore caution that NED redshifts 
should be considered the least reliable of our counterpart as¬ 
signments. The NED associations are dominated by optical 
associations with Abell (Abell 1958) and Northern Optical 
Cluster Survey clusters (Gal et al. 2003). 

Einally, we add four high-z counterparts confirmed using 
SDSS data and which are discussed in Appendix A. 

® The NASA/IPAC Extragalactic Database (NED) is oper¬ 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 
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7.4. IR information 

At Planck detection positions, we have searched for galaxy 
overdensities in the AllWISE mid-infrared source catalogue 
(Cutri et al. 2013). The AllWISE source catalogue includes 
the combined cryogenic and NEOWISE (Mainzer et al. 
2011) observations from the Wide-held Infrared Survey 
Explorer mission (WISE; Wright et al. 2010). The data 
cover the entire sky and we used the deepest bandpasses, 
the 3.4 /r (Wl) and 4.6 /r (W2) channels. We predicted 
galaxy (W1-W2) colours from Bruzual & Chariot (2003) 
stellar population models, and searched for galaxy over¬ 
densities of the same colour in successive redshift ranges 
from z=0.3 to z=1.5 (e.g. Papovich et al. 2010; Mei et al. 
2012, 2014; Stanford et al. 2014). At redshift z < 0.3, the 
contrast between red mid-infrared galaxies and the back¬ 
ground is not efficient for galaxy cluster detection, so we 
searched only in the Helds of Planck detections already val¬ 
idated at redshift z > 0.3, and detections not yet conhrmed 
or with unknown redshift. 

We estimated a signihcance of the over-densities by 
comparing the number of galaxies found in a region of co¬ 
moving diameter of 1 Mpc, with the background galaxies 
found in regions of the same area. To estimate the back¬ 
ground density, we calculated for each redshift range both 
a local background for each candidate, in a region within 
15 arcmin from the Planck detection, and a master back¬ 
ground derived from the estimators for all the Planck Helds. 
A substantial percentage of the Planck detections (~ 37%) 
are affected by artefacts from bright stars in the WISE data, 
which compromise meaningful assessment of the galaxy 
over-densities. The bulk of these are at low galactic lati¬ 
tude (|6| < 20). This means that we do not expect to reach 
a detection completeness of better of 60 — 70%. 

After visual inspection of all detections, we have Hagged 
our detections in regions not affected by bright star arte¬ 
facts with the following classiHcation: 3 - SigniHcant galaxy 
overdensity detected ; 2 - Probable galaxy overdensity; 1 
- Possible galaxy overdensity; 0 - No significant galaxy 
overdensity. We also include these classifications for detec¬ 
tions in regions affected by bright star artefacts: -1 Possible 
galaxy overdensity; -2 - No significant galaxy overdensity; 
-3 No assessment possible. 

To test our classification, and evaluate our completeness 
and purity, we blindly apply our automated and visual in¬ 
spection to 100 objects: 50 confirmed z > 0.5 clusters and 
50 random positions in the sky. We show the results of this 
validation test in Table 6: 59% of the fields have images 
with bright star artefacts, including 17% which are class 
-3. In class 3 and 2, we obtain a 96% and 80% purity, re¬ 
spectively, for the validated clusters. Given the high purity 
of the class 3 detections, we classify these objects as con¬ 
firmed infra-red clusters. 

In Table 7, we show the number of WISE Planck de¬ 
tections in each class for the 935 Planck detections with 
z > 0.3 or unknown redshift. A detailed study of the WISE 
detections will be published in a separate paper (Mei et al. 
2015, in prep.). 73 new clusters have been confirmed (class 
3 and not validated by other methods) by our WISE image 
analysis. A further 54 probable new clusters are identified 
(class 2 and unvalidated by other methods). 


Table 6. Total, validated, and spurious detections in a 
blind test of our WISE detection classification using 50 real 
and 50 spurious fields. 


Test Detections 


Class Total Validated Spurious 


3 . 24 23 1 

2 . 10 8 2 

1. 3 0 3 

0. 4 0 4 

-1. 14 2 12 

-2 . 28 8 20 

-3. 17 9 8 


7.5. SZ information 

We searched for counterparts of our Planck detections 
using catalogues obtained with other SZ surveys, such 
as the South Pole Telescope (SPT) and the Atacama 
Cosmology Telescope (ACT), and update the list of SZ 
confirmations by direct follow-up with the Arc-minute 
Micro-kelvin Interferometer (AMI). 


7.5.1. SPT 

For SPT, we refer to the recently released catalogue 
(Bleem et al. 2014), extracted from the full 2500 deg^ SZ 
survey. It contains 677 SZ detections, of which 516 have 
been confirmed as clusters, through optical and near-IR 
observations. The catalogue contains the photometric 
redshifts (spectroscopic when available) and mass estimate 
of the confirmed clusters. More specifically, we chose to 
use the “Fiducial Cosmology” Catalogue provided by the 
SPT collaboration, to be consistent with the cosmological 
parameters used in the present paper. We performed a 
two step matching process as described in Sect. 7.2: a 
positional match within 10 arcmin, leading to 89 single and 
five double matches, which we verified by comparing our 
mass estimate with the one provided in the SPT catalogue. 
The mass estimates are usually consistent at better than 
3cr, except for one detection, PSZ2 G249.87 — 21.65, 
where they differ at 3.5cr. We note, however, that the 
error bars on the Planck mass reflect only the statistical 
error on the probability contours and do not consider the 
uncertainties nor the intrinsic scatter in the scaling relation 
used to break the degeneracy. Moreover, as discussed also 
in the SPT case by Bleem et al. (2014), the use of a 
fixed scaling relation and of a fiducial cosmology results 
in an underestimation of the statistical and systematic 
uncertainties in both datasets. Therefore, we decided to 
keep the SPT counterparts for PSZ2 G249.87 — 21.65, 
which is also associated in the SPT catalogue to the same 
MCXC cluster (RXC J0628.8 — 4143) as in our matching. 
We checked the position of the matches in the separation 
plane: none of the single matches have been discarded in 
this way, while in four out of five multiple matches, one 
of the counterparts was excluded following this criterion. 
In the remaining match, two possible counterparts are 
allowed and we selected as the most likely the one with a 
smaller mass difference and a smaller separation in terms 
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Table 7. WISE Planck detection classification. 


WISE Planck Detections® 


Previously 

Class Total % of Sample Confirmed Unconfirmed 


3 . 374 40 301 73 

2 . 68 7 14 54 

1 . 55 6 7 48 

0 . 88 9 15 73 

-1 . 42 4 5 37 

-2 . 97 10 15 82 

-3 . 211 23 55 156 


“ For each WISE detection class, we show the total number of Planck detections and their percentage with respect to the 935 
objects with known redshift 2 : > 0.3 or unknown redshift, the number of previously confirmed Planck clusters and the number 
of unconfirmed Planck detections. 


of 0500- 

We observe a systematic difference between the masses 
we derived from Planck data and those provided in the 
SPT catalogue. The mass is a derived quantity which re¬ 
quires scaling information to be assumed before it can be 
calculated from the SZ signal measured by either instru¬ 
ment. Comparison of the SZ observables is complicated by 
the different scales probed by each instrument: Planck is 
sensitive to the cluster outskirts while SPT is sensitive to 
the core regions. Any comparison necessarily requires model 
extrapolation, which is complicated further by the different 
pressure models used in the two measurements. A robust 
comparison will require a joint analysis of the data, which 
is beyond the scope of this paper. 

7.5.2. ACT 

For ACT, we use the catalogue published in Hasselfield 
et al. (2013), which contains both the most recent ACT 
detections and an update of the 23 Marriage et al. (2011) 
detections. 32 PSZ2 detections match with ACT clusters in 
a 10 arcmin radius. 28 have ACT UPP-based SZ masses 
consistent with Planck Mgz (less than Scr deviation). 
Four have more than Scr deviation: PSZ2 G053.44-36.25 
(ACT-CLJ2135.1-0102) and PSZ2 G130.21-62.60 (ACT- 
CLJ0104.8-1-0002) are actually considered as good matches 
by our redMaPPer association, PSZ2 G262.27-35.38 (ACT- 
CLJ0516-5430) is also a good match for our MCXC associ¬ 
ation. The last one, PSZ2 G265.86-19.93 (AGT-GLJ0707- 
5522), is a good match in the PSZl. We decided to leave 
blank the AGT field of the PSZ2 catalogue for these four 
clusters but, given the uncertainties on the mass determi¬ 
nation and associated errors, we did not break the corre¬ 
sponding redMaPPer, MCXC and PSZl association. 

7.5.3. AMI 

Following the ongoing follow-up observations of the Planck 
cluster candidates, a total of 161 clusters with 4.5 <S/N< 
20 were observed with AMI. The detection significance is 
then characterised by calculating of the natural logarithm 
of the Bayes factor. Initio, 

Initio = AlnZio = InZi — InZo, (9) 


Table 8. AMI scale for an interpretation of the detection 
significance of the Planck cluster candidates. 


Category 

In Bio 

N 

Clear detection. 

InBio > 3 

102 

Moderate detection. 

0 < InBio < 3 

30 

Non-detection. 

—3 < InBio < 0 

25 

Clear non-detection. 

ln Bio —3 

4 


where ln2^i and ln2^o are the natural logarithm of the 
Bayesian evidence for model Pfi and Pto respectively. Model 
Hi accounts both for the cluster signal and the contribu¬ 
tion from radio sources while Hq only takes into account the 
radio source environment. Further details of the AMI obser¬ 
vations, of the Bayesian methodology and of the modelling 
of interferometric SZ data, primordial GMB anisotropies, 
and resolved and unresolved radio point sources, as well 
as of the criteria used to categorise clusters are given in 
Perrott et al. (2014) and references therein. 

In this context the detection significance for the 161 
clusters is described in Table 8. 132 of the 161 AMI- 
confirmed Planck clusters are included in the PSZ2. 

7.6. Redshift compilation 

We provide a single redshift estimate for each detection 
where at least one redshift is known for a matched coun¬ 
terpart. In many cases, we have multiple estimates per de¬ 
tection. This section discusses the compilation of redshift 
information and how redshifts are assigned in the final cat¬ 
alogue. Gonfirmation statistics and final assigned redshift 
numbers are summarised in Table 9. 

Initial redshift estimates are taken from the PSZl 
catalogue redshift compilation, given the matching in 
Sect. 7.1. We include new follow-up results from the Planck 
collaboration (Planck Gollaboration et al. 2014; Planck 
Gollaboration Int. XXXV 2015), which include spectro¬ 
scopic updates to PSZl photometric redshifts and new con¬ 
firmations of PSZl detections. We also include external 
updates to PSZl counterpart redshifts from the NED and 
SIMBAD databases (Planck Collaboration XXXVI 2015). 

After these steps, we cycle through priority levels in our 
systematic counterpart searches, in the following order of 
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priority: MCXC, redMaPPer ACT and SPT. We compare 
the updated PSZl and MCXC redshifts with the redshifts 
from redMaPPer where available, and prioritise redMaPPer 
redshifts highest amongst available photometric redshifts. 
We test spectroscopic redshifts at z > 0.1 for consistency 
within Az < 0.03 and Az/z < 0.1 of the redMaPPer photo- 
z. Any discrepancies are considered on an individual basis. 
In a small number of cases, we choose the redMaPPer red- 
shift. We also reject a small number of counterpart assign¬ 
ments where that counterpart is a bad match in redMaPPer. 
These cases are discussed in Appendix B. 

The common sample between the PSZ2 and each of 
the external samples is denoted in the catalogue. After the 
systematic searches, we assign any remaining unconfirmed 
clusters to database counterparts where available. If a PSZl 
match assigns to the same counterpart as one of the negated 
counterparts from the systematic searches, where the pos¬ 
sible counterpart violated the consistency criteria, then the 
PSZl match is also negated. 

8. Sample Properties 

8.1. Mass and redshift properties 

We discuss here the distribution of Planck SZ-selected clus¬ 
ters in the mass-redshift (M 500 — z) plane, using the mass 
proxy derived with scaling relations as discussed in Sect. 7. 
For 1094 detections with known redshifts in the PSZ2 cat¬ 
alogue (Sect. 7), we show in Fig. 26 their position in the 
-^500 ~ plane, compared with the expected completeness 
function C{M^qq,z) of our survey (we show the 20%, 50% 
and 80% completeness levels). These curves indicate the 
points in the M 500 — 2 : plane at which clusters have C% 
chances to be detected. They were computed for the full 
survey area. The red points in Fig. 26 show the 298 new 
PSZ2 confirmed detections, with redshifts, that were not 
found in the previous version of the catalogue. The black 
points show the common PSZ1-PSZ2 detections. 

We stress that the M 500 — 2 distribution in Fig. 26 can¬ 
not be considered as fully representative of the Planck SZ 
selection, since it reflects the biases due to the non-uniform 
knowledge of redshifts over the sky in the ancillary informa¬ 
tion we used (Sect. 7). For instance, we have an extensive 
redshift information in the sky area covered by the SDSS 
survey thanks to the redMaPPer catalogue (Rykoff et al. 

2014) but not in the remaining part of the sky. The incom¬ 
plete redshift information can also explain the rarity of new 
detections in the PSZ2 catalogue with respect to PSZl at 
high redshift: at z > 0.6 we have 36 objects, but only four 
new PSZ2 detections. We note however that most of the 
PSZl clusters in this redshift range were not present in ex¬ 
isting catalogues but they were confirmed as clusters and 
their redshift was measured thanks to the massive follow up 
campaign with optical and X-ray telescopes which was un¬ 
dertaken by the Planck Collaboration for PSZl candidates 
(Planck Collaboration XXIX 2014 and references therein) 
and which continued also after the 2013 release (Planck 
Collaboration et al. 2014; Planck Collaboration Int. XXXV 

2015) . Since a similar observational campaign has not yet 
been possible for new PSZ2 detections, we could not pop¬ 
ulate further the high-mass high-z part of the M 500 — 2 
plane. 

The new PSZ2 confirmed detections (red points in 
Fig. 26) are mostly low-mass objects close to the detection 



Fig. 27. Box-and-whisker diagrams showing the mass dis¬ 
tribution of the PSZ2 (red), PSZl (blue) and their intersec¬ 
tion (black) sample in seven redshift bins. The bottom and 
the top of the boxes represent the first and third quantile of 
the data, while the band inside the box shows the median 
(i.e. the second quantile). The ends of the whiskers mark 
the minimum and maximum of the data. 


limit of the survey. The mean mass of confirmed clusters 
over the whole redshift range in the PSZ2 is 4.82 x 10^^ Mq, 
which is lower than in the PSZl (5.12 x 10^"^ Mq). The com¬ 
mon sample of 795 objects contains the higher mass clusters 
detected by both surveys, with mean mass 5.16 x 10^"'^ Mq. 
This is expected, since the common sample contains none 
of the new low-mass PSZ2 detections and none of the miss¬ 
ing low S/N PSZl detections, discussed in Sect. 6.3, which 
were likely to have been low mass. 

This is also shown in Fig. 27, where we compare the 
mass distribution of the confirmed clusters in the PSZ2, 
the PSZl and their common sample, for several redshift 
bins. The median mass and the first and third quantiles 
are always lower for the PSZ2 than for the PSZl and the 
common sample, showing that we are significantly expand¬ 
ing the sample towards lower masses. 

Fig. 26 also shows a comparison of the SZ selected sam¬ 
ples from the Planck., ACT and SPT surveys. Planck tends 
to detect the rarest high-mass clusters observed at high- 
redshift in these partial-sky surveys and provides a com¬ 
plementary clean SZ selection at lower redshifts, where the 
Planck frequency range provides sufficient information to 
disentangle the SZ signal of large clusters from the back¬ 
ground. 

8.2. X-ray underluminous clusters 

The presence of a bright cool-core, characterised by a 
peaked surface brightness profile, has been shown to bias X- 
ray flux selected cluster samples in favour of peaked, relaxed 
objects with respect to morphologically disturbed systems 
(Eckert et al. 2011). In contrast, SZ selected samples have 
produced more disturbed systems than expected, with SZ 
discovered clusters typically lying on the lower end of the 
mass-luminosity relation (Planck Collaboration IX 2011). 
There has also been much interest in the possible existence 
of severely X-ray under-luminous clusters, with several au¬ 
thors identifying potential systems (eg: Bower et al. 1997; 
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Table 9. Summary of ancillary information. The highest available priority redshift source, following the ordering in 
the Priority column, provides the reference confirmation and redshift. When two priorities are given, the first number 
pertains to spectroscopic redshifts and the second number to photometric redshifts. The PSZ2 contains 1203 confirmed 
clusters, of which 289 are Planck-discovered. 87 of these are clusters newly identified in this paper: 73 are confirmed by 
WISE, eight are new identifications in SDSS data and six are confirmed by AMI. 


Confirmation 

source 

Validation 

Priority 

Joint 

sample size 

Reference 
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Fig. 26. Left panel: Distribution of the 1094 PSZ2 clusters with counterparts with known redshift in the M 500 — z plane. 
New PSZ2 detected clusters are indicated with red dots, while commmon PSZl and PSZ2 clusters are indicated by black 
dots. The solid, dashed and dotted lines indicate respectively the 20%,50% and 80% survey completeness contours for 
the PSZ2. Right panel: Distribution of the PSZ2 clusters with associated redshift in the M 500 — z plane compared to the 
SPT (Bleem et al 2014) and ACT (Hasselfield et al 2013) catalogues. Black circles represent PSZ2 clusters, while red and 
green filled circles mark common SPT/PSZ2ACT/ PSZ2 clusters, respectively. The remaining SPT and ACT clusters 
not detected by Planck are shown with red and green empty squares. 


Popesso et al. 2007; Dietrich et al. 2009; Trejo-Alonso et al. 
2014 and references therein) and suggesting a model where 
these clusters are dynamically young objects, still undergo¬ 
ing accretion and mergers and yet to reach equilibrium. An 
alternative suggestion is that line-of-sight structures may 
bias mass and richness estimates high relative to the X-ray 
luminosity (Bower et al. 1997; Giles et al. 2015). 

However, the reported under-luminosity of these objects 
is disputed. Andreon & Moretti (2011) note that the under¬ 
luminosity is often claimed relative to biased scaling rela¬ 
tions, and that the significance of the under-luminosity is 
amplified due to underestimation of the true scatter in the 
relation. 


In the SDSS area, the majority of Planck detections 
possess counterparts in the redMaPPer catalogue, with red¬ 
shifts and optical richness estimates. We construct a test 
sample from the Planck-redMaPPer intersection at low red¬ 
shift, z < 0.2. This sample of 148 clusters can be expected 
to be detectable in the ROSAT maps. While redMaPPer is 
not complete for the Planck mass ranges at these redshifts, 
it allows us to construct a sample that is independent of any 
X-ray selection effects and therefore well-suited for finding 
under-luminous clusters, as any biases in selection will af¬ 
fect both normal and under-luminous X-ray cluster alike. 

For each of these clusters, we calculate the X-ray count- 
rate using a growth-curve analysis on the ROSAT 0.5 — 2.4 
keV band maps following Bohringer et al. (2000). We de- 
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rived the count-rate and its upper and lower limit from the 
growth curve at 05 oo, which we calculated from the Planck 
mass proxy. We then converted the RASS count rates into 
flux in the 0.1 — 2.4 keV energy range, using an absorbed 
thermal model, where we used the galactic absorption at 
the Planck position, the redshift of the redMaPPer coun¬ 
terpart and the temperature derived from the mass proxy 
through the M-T scaling relation by Arnaud et al. (2007). 
We then converted fluxes into luminosity in the 0.1 — 2.4 
keV channel. We found good agreement with the reported 
MCXC values of L 500 for those clusters also present in the 
MCXC. 

We then searched for outlier clusters from the k 5 oo --^500 
and Rsoo—A relations, using the Y 500 calculated following 
Sect. 5.3. In both cases, we find the best-fit relation for 
our sample using the BCES algorithm. To exclude outliers 
from the fits, we clipped objects with orthogonal residual 
|r_L| > 2.5(Ttot from the best-fit relations, where atot is the 
raw scatter around the relation derived from the median- 
absolute deviation. We then iterated the sigma-clipping 
process until converged. 

The top panels of Fig. 28 show the best-fit relations and 
their zL2atot scatters between Ysoo and redMaPPer richness 
A, and between Ysoo and T 500 , for the test sample of 148 
Plancfc-redMaPPer clusters. The bottom panels of Fig. 28 
show the histograms of A, normalised by Utot- 

The points highlighted with larger cyan circles denote 
clusters with T 500 more than 2.5crtot below the best-fit rela¬ 
tion. These clusters are under-luminous in X-rays for their 
Fsoo- However, their Y 500 estimates are consistent with their 
optical richness and do not lie preferentially beneath the re¬ 
lation. The consistency of optical and SZ mass proxies sug¬ 
gest either that these clusters are under-luminous for their 
mass, or that both the Tsoo and A estimates are biased high. 

One exception to this is PSZ2 G127.71-69.55, discrepant 
with both relations: at lOctot for ^ 500—^500 and Satot 
for Ysoo—A. It is circled in red in Figure 28. This is the 
only cluster in the Planck-redKaPPer sample with a poor 
Q_NEURAL flag and is either a failed redMaPPer match 
(lying close to the matching threshold in the Msz — A plane) 
or a cluster with a severely IR contaminated spectrum for 
which the Planck Ysoo estimate is likely overestimated. We 
therefore remove this from the list of under-luminous can¬ 
didates. 

We expanded the search for under-luminous clusters 
across the whole sky, testing all clusters with z < 0.2 
against the mean hsoo—T 500 relation determined above. 
This expanded the number of under-luminous candidates 
to 22. These objects warrant follow-up observations to de¬ 
termine their dynamical state and to search for line-of-sight 
structures that may bias high both the T 500 and optical 
richness. They are flagged in the COMMENT field of the 
catalogue. They should be of interest in understanding the 
systematic differences between SZ and X-ray selection. 

9. Summary 

The Planck satellite is unique in providing broad frequency 
coverage over the whole sky with good sensitivity to both 
the high frequency spectral increment and the low fre¬ 
quency decrement of the thermal SZ effect. In this paper, 
we have presented second Planck catalogue of SZ sources 
(PSZ2). This is based on data from the full 29 month mis¬ 
sion and uses a methodology that refines the one used to 


produce the PSZl from 15.5 months of data. The catalogue 
is based on the union of results from three cluster detection 
codes (Planck Collaboration XXIX 2014). The PSZ2 con¬ 
tains 1653 cluster candidates distributed across 83.6% of 
the sky. The catalogue was validated using external X-ray, 
optical, SZ and near infra-red data, producing confirmation 
for 1203 candidates with 1094 redshifts. The catalogue con¬ 
tains 716 new detections including 366 confirmed clusters 
with newly identified SZ signal. 87 of our confirmed clusters 
are newly identified in this paper. We have found good con¬ 
sistency with the PSZl and re-detect 937 SZ sources from 
the PSZl sample of 1227. We have investigated the missed 
detections: the vast majority of these were low-significance 
PSZl detections whose S/N has fluctuated beneath the de¬ 
tection threshold. The majority of these are expected to be 
spurious detections. 

The current status of our knowledge of counterparts 
for our detections at various frequencies is summarised in 
Table 10 and compared to the PSZl. Our optical valida¬ 
tion scheme is based on the newly released SDSS-based 
redMaPPer catalogue (Rozo & Rykoff 2014). This produced 
374 high-quality matches where the counterparts are consis¬ 
tent with Planck mass and position information. We reject 
188 possible matches where the mass or position informa¬ 
tion is inconsistent with the Planck information. This un¬ 
derlines the importance of consistency checks when match¬ 
ing with high density SDSS catalogues. Our X-ray and SZ 
counterpart searches implement similar consistency criteria 
leading to tight control over mismatches. 

Central to the counterpart search process is the under¬ 
standing of the Planck SZ parameter estimates. We have 
validated the Planck Compton-Y posteriors using detailed 
simulations that include an ensemble of hydro-dynamically 
simulated pressure profiles that vary from the pressure pro¬ 
file assumed by our extraction algorithms. Our Y 5 R 500 es¬ 
timates are robust to mis-matches in the pressure profile. 
When translating to T 500 , we have shown the importance 
of accurate prior information about radius to break the 
Y — 9 degeneracy and produce accurate and precise es¬ 
timates from the Planck data. Our counterpart searches 
make extensive use of the Planck mass-proxy; this uses 
prior information about redshift and scaling relations to de¬ 
rive mass constraints which show low scatter with respect to 
external estimates. We provide this Msz for all candidates 
with a redshift, and provide Msz{z) in the range 0 < z < 1 
for all other candidates. We expect this information to be 
useful in future comparisons with external data. 

Central to any statistical use of a cluster sample is the 
survey selection function. We have estimated the catalogue 
completeness using Monte-Carlo source injection and we 
provide this as a product for the full survey and for various 
sub-samples as a function of selection S/N. We have vali¬ 
dated the completeness through a comparison with external 
X-ray data and high resolution SZ data from SPT (Bleem 
et al. 2014), which spans the redshift range and angular 
sizes of the Planck data. We estimate the catalogue to be 
83-87% pure, based on simulations of the Planck data and 
detection-by-detection quality assessment utilising machine 
learning. Higher reliability sub-samples can be constructed 
easily: the main contaminant is infra-red galactic emission 
and as such the reliability is a strong function of galactic 
latitude. Specifically, the cluster cosmology zone that covers 
65% of the sky contains 1308 detections at ~ 90% reliabil¬ 
ity, and full survey intersection catalogue (objects detected 
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Fig. 28. Properties in the Psoo—A (top panels) and Psoo-^soo (bottom panels) planes for the Planck-redHaPPer sample 
at z < 0.2. Under-luminous candidates are denoted with cyan circles in the scatter plots to the left, which also show the 
best fit relation and the dispersion ±2tTtot- The circled red point is a cluster with contaminated Y signal. The right plots 
show the histograms of orthogonal deviation Aj_ for each relation. 


by all three codes) contains 827 detections at > 95% relia¬ 
bility. 

Cosmology using the cluster counts is also dependent 
on external observational data to provide cluster redshifts. 
Planck Collaboration XXIV (2015) have produced cosmo¬ 
logical constraints using samples drawn from the PSZ2, 
containing 493 candidates from the intersection sample 
and 439 drawn from the single-code MMF3 sample. Utilising 
larger samples from the PSZ2 requires further redshift in¬ 
formation. We also expect the PSZ2 to contain many high- 
mass clusters at z > 0.6. So far only 36 have been identi- 
hed, of which 21 were identified in targeted follow-up ob¬ 
servations of PSZl candidates. For these reasons, the PSZ2 
should motivate further follow-up observations. In particu¬ 
lar, the catalogue contains 73 clusters confirmed by WISE 
infra-red data that currently have no redshift information 
but which are likely to be high-redshift. 

Understanding the biases in cluster selection that af¬ 
fect samples defined at different wavelengths will be impor¬ 
tant for interpreting statistical results from existing surveys 
and those planned for the near future. Using a low-redshift 
overlap sample from PSZ2 and redMaPPer, we have iden- 


tihed a population of low-z clusters with ‘typical’ optical 
and SZ properties, but which are underluminous for their 
mass in ROSAT X-ray data. These clusters may be part 
of a population of dynamically disturbed clusters that are 
under-represented in X-ray selected surveys. These objects 
will be interesting targets for multi-wavelength follow-up to 
determine their dynamical state. 

In the near future, Planck all-sky SZ data can be com¬ 
bined with observations of the large-scale structure by sur¬ 
veys such as PAN-STARRS, LOFAR, Euclid, LSST, and 
RSG/e-ROSITA. This will provide an unprecedented multi¬ 
wavelength view of the evolution of large-scale structure 
that will revolutionise our understanding of the physics gov¬ 
erning this process. 
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Table 10. Counterpart summary for PSZ2 compared to PSZl. Common samples are defined as those PSZ2 detections 
with the given property that has a counterpart with that property in the PSZl 2015. The intersection comprises those 
detections common to all three detector codes. Low reliability candidates possess a poor neural-network quality assessment 
flag. In the PSZl, low reliability candidates possess the lowest external quality assessment flag. SZ clusters denote clusters 
with SZ detections in ACT or SPT. PSZl 2013 refers to the 2013 release of the catalogue (Planck Collaboration XXIX 
2014), and PSZl 2015 to a recent addendum updating the counterpart information of the catalogue (Planck Collaboration 
XXXVI 2015). 


Sample 

PSZl 2013 

PSZl 2015 

PSZ2 

Common 

New PSZ2 

Union. 

1227 

1227 

1653 

937 

716 

Intersection. 

546 

546 

827 

502 

325 

Confirmed. 

861 

947 

1203 

820 

383 

Candidates . 

366 

292 

546 

99 

447 

Low reliability .... 

142 

131 

143 

39 

104 

Total X-ray. 

501 

501 

603 

477 

126 

MCXC. 

455 

455 

551 
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SZ clusters. 

82 

82 

110 

79 
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Appendix A: High-redshift SDSS confirmations 

Table A.l gives optical information for four high-redshift 
confirmations found using a search in SDSS data around 
unmatched Planck detections. 

We use a multi-wavelength approach to confirm the 
clusters. Each of these candidates possess coincident high- 
redshift optical over-densities in SDSS, firm infra-red confir¬ 
mations from WISE, and significant emission in the ROSAT 
0.5 — 2.4 keV band. We estimate the X-ray luminosity from 
the ROSAT maps using growth-curve analysis, and confirm 
that the luminosity is consistent with the measured Mgz as 
discussed in Sect. 7.2. 

One interesting case is PSZ2 G097.52-1-51.70, which ap¬ 
pears to be a near line-of-sight projection with components 
at z = 0.7 and z = 0.333, separated by 1.91 arcmin. Both 
systems may contribute to the observed Planck signal. We 
have associated to the high-redshift cluster because it is 
significantly closer to the SZ centre (0.714 arcmin vs. 2.47 
arcmin separation), and because it is coincident (0.23 ar¬ 
cmin) with the ROSAT X-ray centre. The z = 0.333 system 
shows less significant hard-band emission. 
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Table A.l. Optical information for our high-z SDSS confir¬ 
mations. Alongside the redshift, we give the RA and DEC 
of the BCG and if the redshift is spectroscopic, gives 
number of cluster members with spectroscopic redshifts. 


NAME 

OIBCG 

Sbcg 

z 

N. 

PSZ2 G076.18-47.30 

343.1475 

4.5381 

0.666 

3 

PSZ2 G087.39-P50.92 

231.6383 

54.1520 

0.748 

1 

PSZ2 G089.39-P69.36 

208.4382 

43.4843 

0.68 


PSZ2 G097.52-P51.70 

223.8374 

58.8707 

0.7 



Appendix B: Differences in PSZl and PSZ2 
redshift assignments 

The PSZ2 contains 782 clusters which had redshift esti¬ 
mates in the PSZl. We assign the same redshift in all but 
43 of these cases. In 25 of these cases, there is no significant 
difference, defined by |Az| > 0.03 or |Az|/z > 0.1, between 
the estimates and we have updated PSZl photometric red¬ 
shift from various sources to new estimates from redMaPPer 
or Planck follow-ups. 

We have updated a further seven PSZl photometric red- 
shifts with recent Planck ENO follow-up redshifts (Planck 
Collaboration Int. XXXV 2015) where the redshift has sig¬ 
nificantly changed. Of these, one was from the PSZl SDSS 
search, two were from PanSTARRs and four were from ear¬ 
lier Planck photometric follow-ups. These updates are in¬ 
cluded in the 2015 update of the PSZl redshift compilation 
(Planck Collaboration XXXVI 2015). 

We discuss the remaining II significant differences be¬ 
low: 

- PSZ2 G020. 66+37.99 (PSZl INDEX 51): 

Rozo et al. (2014) discuss this cluster in depth. 
RedMAPPer finds two overlapping clusters, Zgpec = 
0.338, A = 85.4 and Zspec = 0.443, A = 23.5. The PSZl 
redshift +hot = 0.39 from WEI Planck follow-up is 
likely to be biased by members from the less rich and 
more distant system. For the PSZ2, we choose the higher 
richness match and quote the Zphot = 0.345. The 2015 
update to the PSZl adopts this change. 

- PSZ2 G066.68+68.44 (PSZl INDEX 222): 

This is an ambiguous system. The PSZl used ^spec = 
0.1813 from NORAS (Bohringer et al. 2000). For the 
PSZ2 we have quoted Zphot = 0.163 from redMaPPer 
which agrees with the estimate Zspec = 0.16 from BCS 
follow-up (Struble & Rood 1999; Crawford et al. 1995) 
and the SDSS BCG estimate Zgpec = 0.163. The decision 
between these two ROSAT follow-up spectroscopic red¬ 
shifts rests on the redMaPPer information, which iden¬ 
tifies a rich A = 84.1 system at Zphot = 0.163. 

- PSZ2 G087.39+50.92 (PSZl INDEX 299): 

The PSZ2 position has moved closer to a clear high- 
redshift SDSS cluster at 2:spec = 0.748 at separa¬ 
tion 0.9 X 0err, and away from the PSZl SDSS match 
which is now at 2.62 x 0err- There is also clear ROSAT 
0.5 — 2.4 keV X-ray emission at the high-z location, 
whose strength is consistent with the SZ emission, while 
there is no significant emission at the PSZl SDSS match 
location. 

- PSZ2 G090.66-52.34 (PSZl INDEX 308): 

The PSZl redshift Zgpec = 0.1784 came from a single 
galaxy spectrum (Struble & Rood 1999). redMaPPer 


suggest this galaxy is likely to be in the foreground, 
with the rich A = 85 cluster at slightly high redshift 
-Zphot = 0.197. We note however that the difference is 
small (10.2%). 

- PSZ2 G113.91-37.01 (PSZl INDEX 416): 

We adopt the Rozo et al. (2014) update of the NORAS 
redshift, which replaced a A = 7.1 group at z = 0.135 
with a rich A = 159 cluster at z = 0.371 separated by 8 
arcmin. 

- PSZ2 G121.13+49.64 (PSZl INDEX 443): 

We note that this system is a probable projection. We 
adopt the correction of the redshift from (Rozo et al. 
2014), noting that the richness of the z = 0.22 compo¬ 
nent is consistent with the SZ signal, while the z = 0.438 
system matched in the PSZl is insufficiently rich. The 
2015 update to the PSZl adopts this change. 

- PSZ2 G143.26+65.24 (PSZl INDEX 513): 

The PSZl association with AGO 1430 is correct, but 
we update the redshift of Zgpec = 0.211 from two mem¬ 
bers (Struble & Rood 1991) with Zphot = 0.363 from 
redMaPPer. The X-ray and optical images show an E- 
W elongation and two possible galaxy concentrations, 
possibly a projection. The high-redshift component has 
richness consistent with the SZ signal. The 2015 update 
to the PSZl adopts this change. 

- PSZ2 G151.19+48.27 (PSZl INDEX 537): 

The PSZl association with A0959 is correct. NED lists 
two literature redshifts for this cluster: 0.289 (which we 
adopt) and 0.353 (adopted in the PSZl). This object is 
bimodal in the optical and in the X-ray and is almost 
certainly a projection. redMaPPer suggests association 
with 0.289 component based on consistency of richness 
with the SZ signal. The 2015 update to the PSZl adopts 
this change. 

- PSZ2 G259.30+84.41 (PSZl INDEX 888): 

We adopt the correction of Rozo et al. (2014) of the 
PSZl redshift Zphot = 0.4125 from NSCS, instead 
matching to a clear and rich SDSS cluster within 1 ar¬ 
cmin of the Planck position at Zphot = 0.323. The 2015 
update to the PSZl adopts this change. 

- PSZ2 G310.81+83.91 (PSZl INDEX 1093): 

The PSZ2 matched to an SDSS cluster at Zphot = 0.446. 
redMaPPer finds this potential counterpart to be insuf¬ 
ficiently rich to be detected by Planck at this redshift 
and there is no X-ray emission in ROSAT. The Planck 
detection may be spurious: there are point sources de¬ 
tected at 353 and 545GHz and the neural network qual¬ 
ity assessment suggests a contaminated spectrum. We 
therefore break the association and leave the detection 
unconfirmed. 

- PSZ2 G318.62+58.55 (PSZl INDEX 1123): 

We adopt the correction of Rozo et al. (2014) of the 
PSZl redshift Zgpec = 0.1144. This associates to a differ¬ 
ent cluster at 4.56 arcmin separation, with redMaPPer 
redshifts Zphot = 0.22 (Zspec = 0.233), and richness 
A = 69.3. 

Appendix C: Modifications of the extraction 
algorithms since the PSZl release 

C.l. MMFl 

The MMFl code used for the PSZ2 is the same as for the 

PSZl with the following changes: 
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— Positions estimates are now calculated with sub-pixel 
positioning using posteriors marginalised over all other 
parameters, rather than taking the pixel centre closest 
to the peak. 

— Position error radius estimation has been debugged. 

— The 2D contour grids for (Y, 9) expand dynamically 
if the 91% confidence region is not entirely contained 
within the grid. This expansion is not applied to S/N< 5 
detections. 

C.2. MMF3 

For the PSZ2 release, we made three improvements on our 

MMF3 code: 

— Bright clusters impact the estimation of background. 
For the PSZl we have adopted two different estimators 
of the MMF3 background depending on the S/N of the 
detections. If the S/N was below 20, the theoretical cal¬ 
culation was used (see Eq. 7 in Melin et al. (2006b)) 
while, for S/N greater than 20, we used the standard 
deviation of the filtered map. When the cluster signal is 
subdominant in the map the two estimators return the 
same result but they differ if the cluster is bright (typ¬ 
ically with S/N above 20) due to the contamination of 
the background by the cluster itself. The choice of using 
two estimators has been made to make the MMF3 back¬ 
ground estimate compatible with PwS and MMFl for the 
PSZl. We tested it against the QA after the 2013 release 
and found that it biases significantly MMF3 two dimen¬ 
sional ( 0 s,F 5 R 5 oo) contours with respect to injected clus¬ 
ter size and flux at high S/N. We thus decided to come 
back to the theoretical calculation of the background 
across all the S/N range as in the earlier version of MMF3 
used for the ESZ. This choice fixes the issue with high 
S/N cluster contours in the QA. But it increases the 
S/N value and shifts the (0 s,F5 R5 oo) contours for MMF3 
detections in the PSZ2 release with respect to PSZl for 
the detections with S/N above 20 more than the increase 
expected from the additional integration time. For the 
PSZ2 PwS now estimate the cross-channel covariance 
matrix under the ^native' prescription (see PwS section 
in this appendix). This improvement makes the back¬ 
ground estimate compatible with the theoretical calcu¬ 
lation from MMF3 and also gives unbiased estimates for 
the PwS two dimensional in the QA. 

— MMF3 two dimensional (^sjbsRsoo) contours in the PSZl 
were tested against the QA after the 2013 release. This 
could not be done in early 2013 for lack of time. The 
contours were found to be wider than expected. The 
code has been corrected and new contours have been 
produced and included in the Planck Legacy Archive 
for the PSZl. The PSZ2 relies on this new and fully 
tested estimate of the two dimensional contours. 

— The MMF3 positional error for the PSZl was overesti¬ 
mated. The code has been corrected and tested against 
the QA. New estimates have been produced and in¬ 
cluded in the Planck Legacy Archive for the PSZl. The 
PSZ2 uses the new estimate of the positional error. 


— The cross-channel covariance matrix is now always esti¬ 
mated using iterative recalibration but is parameterised 
in order to produce a smoother estimate. Using the QA, 
we have shown that the new calibration only impacts 
the S/N estimate (^ 12%) keeping unchanged all other 
parameter estimates. The new S/N estimates are con¬ 
sistent with the other codes. 

— For (0 s,F5 R5 oo)i we now adopt non-informative pri¬ 
ors, formulated using Jeffrey’s method (Carvalho et al. 
2012 ), instead of informative priors derived from a fidu¬ 
cial cosmology and mass function. 


C.3. PwS 

The PwS code used for PSZ2 is similar to the one used for 
PSZl, with two modifications: 
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Appendix D: Description of the delivered products 

The data products comprise: (i) the main catalogue, which contains the characterised catalogue with ancillary information; 
(ii) individual algorithm catalogues produced by each of the codes prior to merging to create the main catalogue, which 
contain Y — 9 parameter posteriors per cluster; (iii) selection function files containing the completeness and survey masks 
for various sample definitions. 

D.l. Main catalogue 

The table contains the following columns: 


Golumn Name 

Data Type 

Units 

Description 

INDEX . 

Integer(4) 


Index of detection (see note 1) 

NAME . 

String 


Name of detection (see note 2) 

GLON . 

Real( 8 ) 

degrees 

Galactic longitude {0° < 1 < 360°) 

GLAT . 

Real ( 8 ) 

degrees 

Galactic latitude (—90° <b< 90°) 

RA . 

Real( 8 ) 

degrees 

Right ascension (J2000) 

DEG. 

Real( 8 ) 

degrees 

Declination (J2000) 

POS.ERR . 

Real(4) 

arcmin 

Uncertainty in position (see note 3) 

SNR . 

Real(4) 


S/N of detection 

PIPELINE . 

Integer(4) 


Pipeline from which information is taken: the reference 
pipeline (see note 4) 

PIPE_DET . 

Integer(4) 


Information on pipelines making detection (see note 4) 

PCCS2 . 

Boolean 


Indicates whether detection matches with any in PCCS2 
catalogues 

PSZ . 

Integer(4) 


Index of matching detection in PSZl, or -1 if a new detec¬ 
tion 

IR FLAG . 

Integer(l) 


Flag denoting heavy IR contamination 

Q_NEURAL . 

Real(4) 


Neural network quality flag (see note 5) 

Y5R500 . 

Real(4) 

10 “^ arcmin^ 

Mean marginal V 5 R 500 as measured by the reference pipeline 

Y5R500_ERR. 

Real(4) 

10 “^ arcmin^ 

Uncertainty onYsRsoo as measured by the reference pipeline 

VALIDATION .... 

Integer(4) 


External validation status (see note 6 ) 

REDSHIFTJD ... 

String 


External identifier of cluster associated with redshift mea¬ 
surement (see note 7). 

REDSHIFT. 

Real(4) 


Redshift of cluster (see note 7) 

MSZ . 

Real(4) 


SZ mass proxy (see note 8 ) 

MSZ_ERR.UP .... 

Real(4) 


Upper la SZ mass proxy confidence interval (see note 8 ) 

MSZ_ERR.LOW .. 

Real(4) 

o 

1—1 

Lower la SZ mass proxy confidence interval (see note 8 ) 

MGXG . 

String 


ID of X-ray counterpart in the MGXG if one is present 

REDMAPPER ... 

String 


ID of optical counterpart in the redMaPPer catalogue if one 
is present 

AGT. 

String 


ID of SZ counterpart in the ACT catalogues if one is present 

SPT . 

String 


ID of SZ counterpart in the SPT catalogues if one is present 

WISE.FLAG . 

Integer(4) 


Confirmation flag of WISE IR overdensity (see note 9) 

AMLEVIDENCE . 

Real(4) 


Bayesian evidence for AMI counterpart detection (see note 

9) ^ 

Indicates whether the cluster is in the cosmology sample 

GOSMO . 

Boolean 


GOMMENT . 

String 


Comments on this detection 


Notes: 

1. The index is determined by the order of the detections in the union catalogue. The matching entries in the individual 
catalogues have the same index to facilitate cross-referencing. 

2. The names are in the format PSZ2 Gxxx.xxiyy.yy where xxx.xx is the Galactic longitude and iyy.yy is the Galactic 
latitude of the detection, both in degrees. The coordinates are truncated towards zero, not rounded. 

3. The value given here is the 95% confidence interval of the distribution of radial displacement. 

4. The PIPELINE column defines the pipeline from which the values in the union catalogue are taken: 1 = MMFl; 2 = 
MMF3; 3 = PwS. 

The PIPEJDET column is used to indicate which pipelines detect this object. The three least significant decimal 
digits are used to represent detection or non-detection by the pipelines. Order of the digits: hundreds = MMFl; tens 
= MMF3; units = PwS. If it is detected then the corresponding digit is set to 1, otherwise it is set to 0. 

5. The neural network quality flag is l-Qbad, following the definitions in Aghanim et al. 2014. Values Q_NEURAL< 0.4 
denote low-reliability detections. 
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6. The VALIDATION column gives a summary of the external validation, encoding the most robust external identifica¬ 
tion: -1 = no known external counterpart; 10=ENO follow-up; 11= RTT follow-up; 12= PanSTARRs; 13= redMaPPer 
non-blind; 14= SDSS high-z; 15=AMI; 16= WISE; 20 = legacy identification from the PSZl 2013 release; 21 = MCXC; 
22= SPT; 23=ACT; 24= redMaPPer; 25= PSZl counterpart with redshift updated in (Planck Collaboration XXXVI 
2015); 30= NED. 

7. The redshift source is the most robust external validation listed in the VALIDATION field. 

8 . Definition of Mgz- The hydrostatic mass, M 500 , assuming the best-fit Y — M scaling relation of Arnaud et al. (2010) 
as a prior that cuts a plane through the parameter contours (see Sect. 5.3). The errors are 67% confidence statistical 
errors and based on the Planck measurement uncertainties only. Not included in the error estimates are the statistical 
errors on the scaling relation, the intrinsic scatter in the relation, or systematic errors in data selection for the scaling 
relation fit. 

9. WISE confirmation flag is assigned by visual inspection and defined to be one of [—10, —2, —1, 0,1, 2, 3], where -10 
denotes no information and the other values are discussed in Section 7.4. Bayesian evidence for AMI counterpart 
defined in paper. 

D.2. Individual algorithm catalogues 

The table contains the following columns: 


Column Name 

Data Type 

Units 

Description 

INDEX. 

Integer (4) 


Index of detection (see note 1) 

NAME . 

String 


Name of detection (see note 1) 

GLON . 

Real(8) 

degrees 

Galactic longitude {0° < 1 < 360°) 

GLAT . 

Real(8) 

degrees 

Galactic latitude (—90° <b< 90°) 

RA . 

Real(8) 

degrees 

Right ascension (J2000) 

DEC . 

Real(8) 

degrees 

Declination (J2000) 

POS.ERR. 

Real(4) 

arcmin 

Uncertainty in position (see note 2) 

SNR. 

Real(4) 


S/N of detection (see note 3) 

TS.MIN . 

Real(4) 

arcmin 

Minimum 9s in second extension HDU (see note 4) 

TS.MAX . 

Real(4) 

arcmin 

Maximum 9s in second extension HDU (see note 4) 

Y.MIN . 

Real(4) 

arcmin^ 

Minimum I 5 R 500 in second extension HDU (see note 4) 

Y.MAX . 

Real(4) 

arcmin^ 

Maximum I 5 R 500 in second extension HDU (see note 4) 


Notes: 

1. The index and name are taken from the union catalogue. The matching entries in the individual catalogues have the 
same index and name to facilitate cross-referencing. 

2. The value given here is the 95% confidence interval of the distribution of radial displacement. 

3. The SNR column contains the native signal-to-noise ratio determined by the detection pipeline. 

4. These entries define the limits of the grid used to evaluate the 2D probability distribution of 9s and V 5 R 500 in the 
second extension HDU (see below). 


Second extension HDU 

The second extension HDU contains a three-dimensional image with the two-dimensional probability distribution in 9s 
and VsRsoo for each cluster. The probability distributions are evaluated on a 256 x 256 linear grid between the limits 
specified in the first extension HDU. The limits are determined independently for each detection. The dimensions of the 
3D image are 256 x 256 x n, where n is the number of detections in the catalogue. The second dimension is 9s and the 
first dimension is Vsrsoo- 


Third extension HDU 

The third extension HDU contains a three-dimensional image with the Msz observable information per cluster as a 
function of assumed redshift. The image dimensions are 100 x 4 x n, where n is the number of detections in the catalogue. 
The first dimension is the assumed redshift. The second dimension has size 4: the first element is the assumed redshift 
value for the Msz fields. The second element is the Msz lower 67% confidence bound, the third element is the Msz 
estimate and the fourth element is the Msz upper 67% confidence bound, all in units of 1O^^M0. These errors are based 
on the Planck measurement uncertainties only. Not included in the error estimates are the statistical errors on the scaling 
relation, the intrinsic scatter in the relation, or systematic errors in data selection for the scaling relation fit. 

D.3. Selection function file format 

The selection function information is stored in FITS files. The filenames of the catalogues are of the form PSZ2- 
selection_Rx.xx.fits, where x.xx is the release number. 


34 


















Planck Collaboration: Planck Legacy SZ 


First extension HDU 

The first extension HDU contains the survey region, denoted by an iVside = 2048 ring-ordered HEALPix map in 
GALACTIC coordinates. Pixels in the survey region have the value 1.0 while areas outside of the survey region have 
value 0.0. 


Second extension HDU 

The second extension HDU contains a three-dimensional image containing the survey completeness probability distribu¬ 
tion for various thresholds. The information is stored in an image of size 30 x 32 x 12. The first dimension is Isoo, the 
second dimension is 05 oo and the third dimension is the signal-to-noise threshold. The units are percent and lie in the 
range 0-100 and denote the detection probability of a cluster lying within the given Tsoo-^^soo bin. 


Third extension HDU 

The second extension HDU contains the Ysoo grid values for the completeness data cube held in the second extension. It 
has length 30 and spans the range 1.12480 x 10“^ — 7.20325 x 10“^ arcmin^ in logarithmic steps. 


Fourth extension HDU 

The fourth extension HDU contains the 05 oo grid values for the completeness data cube held in the second extension. It 
has length 32 and spans the range 0.9416-35.31 arcmin in logarithmic steps. 


Fifth extension HDU 

The fifth extension HDU contains the signal-to-noise threshold grid values for the completeness data cube held in the 
second extension. It has length 12 and contains thresholds at intervals of 0.5 from 4.5 to 10.0. 
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Appendix E: Detail of missing PSZl detections 


Table E.l. Detail of the 291 PSZl detections not present in the PSZ2 catalogue. The TYPE column lists the reason why 
the detection was dropped. TYPE 1 lost detections are low-S/N detections lost due to changes in the noise realisation. 
The S/N_non_blind field contains the non-blind S/N for the Y signal in the full-mission maps at the location and size 
of the PSZl detection and is provided for all TYPE 1 lost detections (whereas the field S/N is for the PSZl). TYPE 2 
are lost behind the new point source mask. TYPE 3 are cut due to IR contamination. TYPE 4 are cut by internal PwS 
consistency criteria. Each of these types are discussed in Sect. 6.3. 


Index 

S/N 

Pipeline 

PlPE DET 

Id ext 

FIedshift 

S/N (non-blind) 

Type 

4. 

6.04 

3 

101 




3 

8 . 

4.92 

1 

100 



0.34 

1 

9 . 

5.76 

1 

100 




3 

13 . 

4.52 

1 

100 

ZwCl 1454.5-1-0656 

0.429 

3.73 

1 

28 . 

4.70 

2 

10 


0.46 

3.54 

1 

30 . 

4.72 

1 

100 



1.83 

1 

32 . 

4.54 

1 

100 



1.65 

1 

34 . 

4.50 

1 

100 



1.42 

1 

38 . 

4.65 

1 

100 



2.48 

1 

40 . 

4.89 

2 

10 



2.90 

1 

41 . 

4.81 

2 

10 



2.90 

1 

43 . 

4.76 

2 

10 



2.69 

1 

52 . 

5.71 

2 

10 


0.39 


4 

58 . 

4.56 

1 

100 



2.06 

1 

59 . 

4.63 

1 

100 



3.27 

1 

60 . 

4.84 

2 

10 

RXC J1917.5-1315 

0.177 

3.11 

1 

61 . 

5.06 

2 

111 


0.650893 

2.91 

1 

62 . 

4.95 

2 

11 

AGO S 1010 

0.28 

4.02 

1 

66 . 

5.20 

2 

10 



3.19 

1 

68 . 

4.67 

3 

1 



3.45 

1 

77 . 

4.58 

1 

100 

RXC J1453.1-I-2153 

0.1186 

2.96 

1 

82 . 

4.98 

2 

10 



2.84 

1 

83 . 

4.96 

2 

10 



3.53 

1 

84 . 

4.84 

2 

10 



3.80 

1 

86 . 

5.23 

1 

100 



2.82 

1 

89 . 

4.69 

3 

101 

WHL J248.764-l-15.4836 

0.4725 

3.27 

1 

90 . 

4.68 

1 

100 



2.11 

1 

97 . 

4.82 

1 

100 

WHL J252.649-l-16.8253 

0.3612 

2.56 

1 

98 . 

4.70 

3 

1 



3.54 

1 

104 . 

4.54 

2 

10 



1.94 

1 

Ill . 

4.57 

2 

10 



1.97 

1 

112 . 

4.61 

1 

100 



1.59 

1 

121 . 

4.70 

2 

no 



1.22 

1 

126 . 

4.67 

2 

10 



1.91 

1 

128 . 

4.54 

3 

1 

RXC J1623.5-I-2634 

0.4274 

3.79 

1 

131 . 

5.23 

2 

11 

AMF J320.551-6.81740 

0.5344 

3.59 

1 

136 . 

4.74 

1 

100 




3 

142 . 

4.52 

3 

101 



1.26 

1 

157 . 

4.64 

1 

100 



2.27 

1 

158 . 

5.29 

2 

10 



2.08 

1 

162 . 

4.60 

1 

100 



2.34 

1 

165 . 

5.11 

3 

101 




3 

170 . 

5.33 

2 

11 



4.32 

1 

175 . 

4.82 

1 

100 


0.1944 

2.84 

1 

176 . 

5.72 

2 

10 




3 

184 . 

4.76 

3 

1 



2.67 

1 

193 . 

4.55 

1 

100 



2.83 

1 

199 . 

4.55 

3 

1 



2.12 

1 

203 . 

4.65 

2 

10 



1.98 

1 

211 . 

4.58 

3 

1 



3.31 

1 

212 . 

4.54 

3 

1 



3.72 

1 

213 . 

4.70 

2 

10 



2.93 

1 

223 . 

4.51 

2 

10 


0.3341 

3.16 

1 

233 . 

4.79 

2 

10 

ZwCl 2151.0-1-1325 

0.205 

3.97 

1 

237 . 

4.78 

1 

100 

AGO 2429 


3.36 

1 

251 . 

5.17 

1 

100 



2.59 

1 

257 . 

4.90 

2 

11 



3.68 

1 

260 . 

4.78 

2 

no 

WHL J242.728-1-51.2267 

0.4096 

3.60 

1 

262 . 

4.52 

1 

100 



3.63 

1 

267 . 

6.37 

2 

111 

ACTJ2327.4-0204 

0.705 


2 
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Table E.2. Continuation of Table E.l. 


Index 

S/N 

Pipeline 

PlPE DET 

Id ext 

Redshift 

S/N (non-blind) 

Type 

271 . 

4.71 

2 

10 



3.84 

1 

272 . 

4.57 

3 

1 

ZwCl 1746.2+5429 

0.31 

4.00 

1 

276 . 

4.65 

2 

11 



2.88 

1 

278 . 

4.55 

2 

10 


0.306807 

2.32 

1 

300 . 

5.07 

2 

10 


0.1132 

4.09 

1 

305 . 

4.82 

2 

10 



2.12 

1 

306 . 

4.90 

1 

100 



3.51 

1 

309 . 

5.37 

2 

110 

ZwCl 1602.3+5917 

0.2544 

4.01 

1 

311 . 

4.64 

2 

110 




3 

314 . 

5.18 

1 

100 



1.78 

1 

317 . 

5.46 

1 

100 



1.09 

1 

321 . 

4.58 

1 

100 

ZwCl 1604.4+6113 

0.3447 

4.06 

1 

327 . 

5.78 

3 

101 

RXC J2318.4+1843 

0.0389 


4 

331 . 

4.79 

3 

101 



2.53 

1 

333 . 

4.71 

3 

1 

WHL J286.905+64.5511 

0.3561 

3.89 

1 

336 . 

6.55 

2 

111 




3 

349 . 

4.87 

2 

11 




2 

361 . 

4.71 

2 

10 



1.28 

1 

365 . 

4.82 

3 

101 

RXC J1834.1+7057 

0.0824 

3.43 

1 

367 . 

4.52 

2 

10 

ZwCl 1748.0+7125 


2.72 

1 

370 . 

5.08 

2 

110 




3 

371 . 

4.55 

2 

110 




3 

372 . 

4.63 

2 

11 



4.03 

1 

373 . 

4.50 

2 

10 




2 

375 . 

4.78 

3 

101 



3.03 

1 

376 . 

4.54 

2 

10 

AMF J359.521+15.1625 

0.1785 

2.56 

1 

382 . 

4.73 

2 

110 




3 

387 . 

4.66 

2 

111 




2 

396 . 

4.55 

2 

10 


0.25 

1.78 

1 

397 . 

6.89 

2 

111 




3 

398 . 

4.55 

2 

10 




3 

400 . 

4.58 

2 

111 


0.533998 

2.74 

1 

405 . 

4.81 

1 

100 

WHL J358.170+38.9803 

0.27 

2.79 

1 

412 . 

4.76 

3 

1 



3.08 

1 

426 . 

5.85 

2 

111 




2 

430 . 

4.54 

3 

1 



3.60 

1 

436 . 

4.96 

3 

101 



3.38 

1 

437 . 

4.51 

1 

100 



3.42 

1 

444 . 

4.72 

2 

110 



3.54 

1 

445 . 

4.66 

2 

10 

Abell 98S 

0.104 

2.87 

1 

446 . 

5.35 

2 

11 




2 

456 . 

4.89 

1 

100 



3.08 

1 

458 . 

4.86 

3 

101 



2.82 

1 

462 . 

5.40 

2 

110 




2 

466 . 

5.91 

2 

11 




3 

468 . 

4.95 

2 

10 



3.00 

1 

469 . 

4.90 

2 

11 


0.423 

2.90 

1 

476 . 

4.91 

2 

10 



2.66 

1 

478 . 

5.42 

3 

1 




3 

479 . 

4.81 

2 

11 



4.25 

1 

483 . 

4.57 

2 

10 



3.40 

1 

488 . 

4.58 

1 

100 

RXC J0115.2+0019 

0.045 

3.60 

1 

489 . 

4.58 

3 

1 



2.75 

1 

490 . 

5.22 

2 

10 



3.18 

1 

491 . 

4.51 

1 

100 

RXC J0152.9+3732 

0.2993 

3.62 

1 

504 . 

4.90 

1 

100 



2.99 

1 

505 . 

4.53 

2 

10 


0.172448 

2.88 

1 

517 . 

4.63 

2 

10 



2.33 

1 

522 . 

6.62 

3 

110 




3 

524 . 

4.96 

3 

1 

RXC J0209.5+1946 

0.0657 

3.57 

1 

527 . 

4.52 

1 

100 


0.38508 

1.54 

1 

529 . 

8.40 

2 

111 




2 

534 . 

4.75 

3 

1 




3 

538 . 

4.70 

1 

100 



3.62 

1 

539 . 

5.59 

2 

111 

ACO 307 


3.52 

1 

544 . 

4.79 

2 

10 



3.57 

1 

555 . 

4.51 

3 

1 



2.81 

1 

556 . 

5.28 

2 

111 


0.532786 

3.54 

1 

557 . 

4.76 

1 

100 



1.97 

1 

559 . 

4.67 

1 

100 



2.48 

1 

562 . 

4.64 

2 

10 

RXC J0157.4-0550 

0.1289 

4.35 

1 

564 . 

6.56 

1 

100 




3 

565 . 

5.22 

2 

111 

RXC J0137.4-1259 

0.2143 

5.16 

1 

576 . 

4.51 

1 

100 



2.81 

1 

586 . 

5.32 

2 

10 



3.04 

1 

590 . 

4.89 

2 

11 



4.83 

1 


37 





















































































Planck Collaboration: Planck Legacy SZ 


Table E.3. Continuation of Table E.l. 


Index 

S/N 

Pipeline 

PlPE DET 

Id ext 

Redshift 

S/N (non-blind) 

Type 

592 . 

4.62 

2 

10 

RXC J0822.1+4705 

0.1303 

3.47 

1 

604 . 

4.68 

1 

100 

RXC J0248.0-0332 

0.1883 

3.19 

1 

605 . 

5.01 

3 

101 



2.68 

1 

607 . 

4.94 

1 

100 

RXC J0956.0+4107 

0.587 

2.96 

1 

611 . 

4.82 

2 

111 




3 

612 . 

5.87 

2 

11 




2 

616 . 

5.09 

2 

10 



2.20 

1 

618 . 

5.43 

2 

11 



3.43 

1 

621 . 

4.78 

2 

11 

RXC J0326.8+0043 

0.45 

2.90 

1 

622 . 

5.43 

2 

10 


0.494731 

1.51 

1 

626 . 

5.90 

2 

111 




3 

629 . 

4.51 

1 

100 



2.07 

1 

639 . 

5.10 

3 

1 



4.27 

1 

642 . 

4.65 

1 

100 

WHL J164.029+34.0043 

0.3805 

3.50 

1 

645 . 

4.51 

1 

100 



2.60 

1 

650 . 

4.56 

1 

100 


0.47 

1.10 

1 

652 . 

4.62 

1 

100 



2.59 

1 

653 . 

5.19 

3 

1 




3 

658 . 

4.57 

1 

100 



1.78 

1 

659 . 

4.82 

1 

100 



3.24 

1 

669 . 

4.85 

2 

10 

RXC JlllO.7+2842 

0.0314 

2.96 

1 

670 . 

4.86 

1 

100 

WHL J56.0261-13.5512 

0.5757 

3.73 

1 

671 . 

4.70 

3 

1 




2 

672 . 

4.94 

2 

110 

WHL J161.821+27.9906 

0.4333 

3.28 

1 

678 . 

4.54 

1 

100 


0.37572 

2.01 

1 

679 . 

4.55 

1 

100 



3.31 

1 

683 . 

5.38 

2 

10 




2 

684 . 

6.62 

3 

101 




2 

697 . 

4.62 

3 

1 

ACO 457 


3.81 

1 

698 . 

4.72 

2 

10 




2 

699 . 

4.77 

1 

100 



2.73 

1 

704 . 

4.56 

1 

100 

WHL J131.956+13.5279 

0.3487 

3.79 

1 

705 . 

4.67 

1 

100 



0.27 

1 

712 . 

4.65 

1 

100 



1.52 

1 

719 . 

4.82 

2 

10 

WHL J158.665+20.5346 

0.4674 

2.95 

1 

721 . 

4.53 

1 

100 

ACO S 270 


2.01 

1 

722 . 

4.71 

3 

1 



2.76 

1 

725 . 

4.56 

2 

10 



3.61 

1 

728 . 

5.29 

1 

100 

RXC J0906.4+1020 

0.1328 

3.67 

1 

729 . 

4.62 

1 

100 

WHL J140.630+11.6581 

0.2609 

2.88 

1 

735 . 

4.52 

3 

1 



3.37 

1 

736 . 

5.02 

2 

10 



3.95 

1 

737 . 

4.52 

2 

10 



3.19 

1 

743 . 

4.57 

2 

10 



3.45 

1 

748 . 

4.55 

3 

1 



3.84 

1 

749 . 

4.97 

3 

101 

ACO 3218 


4.03 

1 

750 . 

5.13 

2 

10 


0.31 

3.37 

1 

751 . 

4.81 

1 

100 

CXOMP J091126.6+055012 

0.7682 

3.82 

1 

753 . 

4.52 

3 

1 



3.42 

1 

755 . 

7.77 

3 

1 




2 

760 . 

4.62 

2 

110 

WHL J134.086+1.78038 

0.7243 

2.90 

1 

762 . 

4.93 

3 

1 




4 

766 . 

7.83 

3 

1 




4 

770 . 

4.64 

3 

1 

RXC J1047.5+1513 

0.2108 

3.50 

1 

771 . 

4.52 

2 

10 

WHL J124.638-6.42296 

0.5123 

3.19 

1 

775 . 

6.97 

3 

1 




2 

781 . 

4.79 

3 

1 




3 

782 . 

4.62 

3 

1 

ZwCl 0919.7-0016 

0.3538 

2.64 

1 

788 . 

4.64 

2 

10 

ACO S 403 


3.72 

1 

789 . 

5.22 

2 

110 




3 

792 . 

4.52 

2 

10 



1.39 

1 

794 . 

5.83 

1 

100 




3 

795 . 

4.60 

2 

111 

WHL J170.480+15.8014 

0.5593 

2.95 

1 

798 . 

5.05 

2 

10 



3.29 

1 

809 . 

4.59 

1 

100 



2.85 

1 

811 . 

4.64 

3 

1 




3 

813 . 

4.64 

1 

100 



3.06 

1 

814 . 

4.58 

3 

1 



4.50 

1 

820 . 

4.57 

1 

100 

RXC J1013.7-0006 

0.0927 

2.49 

1 

827 . 

8.14 

3 

1 




2 

830 . 

4.56 

2 

10 



2.33 

1 

832 . 

4.55 

2 

10 

ACO S 539 


3.04 

1 

833 . 

4.62 

2 

110 



3.38 

1 

836 . 

4.76 

2 

110 

RXC J0345.7-4112 

0.0603 

3.89 

1 

843 . 

4.51 

3 

1 



2.78 

1 


38 
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Table E.4. Continuation of Table E.l. 


Index 

S/N 

Pipeline 

PlPE DET 

Id ext 

Redshift 

S/N (non-blind) 

Type 

844 . 

4.78 

3 

1 



3.91 

1 

845 . 

7.14 

3 

1 




4 

859 . 

4.89 

3 

1 



1.87 

1 

860 . 

4.64 

1 

100 



2.41 

1 

864 . 

5.04 

2 

10 



3.54 

1 

866 . 

5.17 

3 

1 




3 

874 . 

4.55 

1 

100 



3.45 

1 

884 . 

4.67 

1 

100 



0.99 

1 

885 . 

4.53 

3 

1 



3.01 

1 

886 . 

4.72 

2 

111 



2.50 

1 

900 . 

4.60 

3 

1 



2.92 

1 

908 . 

4.68 

2 

10 


0.45 

3.14 

1 

909 . 

4.68 

1 

100 



1.94 

1 

913 . 

4.87 

2 

10 



3.50 

1 

917 . 

4.78 

1 

100 



1.46 

1 

921 . 

4.72 

2 

10 


0.26 

4.23 

1 

923 . 

4.71 

1 

100 



1.57 

1 

925 . 

5.00 

1 

100 



2.16 

1 

927 . 

4.72 

1 

100 



1.84 

1 

928 . 

4.66 

3 

1 



2.65 

1 

933 . 

4.60 

2 

10 



3.06 

1 

949 . 

4.86 

3 

1 



1.88 

1 

950 . 

4.76 

1 

100 



2.11 

1 

953 . 

4.87 

3 

1 



3.98 

1 

964 . 

4.95 

2 

10 


0.14 

2.83 

1 

965 . 

4.51 

1 

100 



1.71 

1 

966 . 

4.51 

3 

1 



2.19 

1 

968 . 

4.54 

2 

10 



2.67 

1 

973 . 

4.56 

1 

100 



1.52 

1 

980 . 

12.78 

2 

111 

RXC J1217.6+0339 

0.0766 


2 

992 . 

5.52 

2 

11 



0.68 

1 

1010. 

5.91 

2 

11 




3 

1016. 

4.76 

2 

110 



2.74 

1 

1018. 

4.61 

3 

1 



3.97 

1 

1019. 

4.98 

2 

11 



3.89 

1 

1031. 

4.81 

2 

110 




3 

1039. 

4.50 

2 

10 



3.30 

1 

1048. 

4.90 

2 

10 

AGO S 137 

0.02764 

2.50 

1 

1049. 

4.95 

2 

10 

AGO 1603 

0.1314 

2.86 

1 

1052. 

5.74 

3 

1 




3 

1055. 

4.62 

3 

1 

RXG J0052.7-8015 

0.1141 

4.03 

1 

1059. 

4.68 

1 

100 



2.74 

1 

1060. 

6.81 

2 

111 




3 

1069 . 

4.84 

2 

10 


0.12941 

2.29 

1 

1080 . 

4.61 

1 

100 



2.40 

1 

1081 . 

4.60 

2 

10 




3 

1091. 

5.21 

2 

11 




3 

1092 . 

4.64 

2 

10 



3.08 

1 

1094 . 

4.68 

1 

100 



2.57 

1 

1103. 

6.15 

1 

100 




2 

1107. 

5.02 

2 

10 



2.39 

1 

1111. 

4.83 

2 

10 



2.74 

1 

1119. 

5.48 

2 

10 



3.53 

1 

1132. 

4.51 

1 

100 



1.43 

1 

1133. 

4.73 

1 

100 


0.25 

2.81 

1 

1135. 

4.78 

3 

1 



3.83 

1 

1144. 

4.90 

3 

1 



4.22 

1 

1152. 

4.63 

3 

1 



3.61 

1 

1155. 

4.78 

2 

111 

SPT-GLJ2148-6116 

0.571 

3.59 

1 

1159. 

4.76 

1 

100 



2.15 

1 

1162 . 

4.63 

1 

100 



2.74 

1 

1170 . 

4.65 

2 

10 



3.63 

1 

1173 . 

4.55 

1 

100 



2.24 

1 

1174 . 

6.61 

2 

111 




2 

1177 . 

5.58 

2 

10 




3 

1178. 

4.50 

1 

100 



2.90 

1 

1180. 

5.88 

2 

111 




3 

1188 . 

4.78 

2 

10 



3.03 

1 

1194 . 

4.75 

2 

10 


0.21 

3.26 

1 

1196. 

4.57 

1 

100 



2.76 

1 

1197. 

4.78 

2 

11 



2.35 

1 

1198. 

4.80 

1 

100 



2.91 

1 

1199. 

4.60 

2 

110 



2.31 

1 

1203. 

5.02 

1 

100 

AGO S 808 

0.049131 

3.41 

1 

1204. 

4.70 

1 

100 


0.5 

3.20 

1 

1212. 

4.74 

1 

100 



2.84 

1 

1215 . 

4.63 

1 

100 



3.43 

1 

1217 . 

4.52 

2 

10 



2.49 

1 

1219 . 

5.16 

3 

1 




2 

1221. 

5.47 

2 

11 




2 
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